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ABSTRACT 


The nearby dwarf starburst galaxy NGC 5253 hosts a number of young, mas¬ 
sive star clusters, the two youngest of which are centrally concentrated and 
surrounded by thermal radio emission (the ‘radio nebula’). To investigate the 
role of these clusters in the starburst energetics, we combine new and archival 
Hubble Space Telescope images of NGC 5253 with wavelength coverage from 
1500 A to 1.9 /im in 13 filters. These include Ha, P/3, and Pa, and the imaging 
from the Hubble Treasury Program LEGUS (Legacy Extragalactic UV Survey). 
The extraordinarily well-sampled spectral energy distributions enable modeling 
with unprecedented accuracy the ages, masses, and extinctions of the 9 opti¬ 
cally brightest clusters (My < —8.8) and the two young radio nebula clusters. 
The clusters have ages rs_/ 1-15 Myr and masses ~ 1 x 10 2 3 4 — 2.5 x 10 5 M 0 . The 
clusters’ spatial location and ages indicate that star formation has become more 
concentrated towards the radio nebula over the last ~ 15 Myr. The most massive 
cluster is in the radio nebula; with a mass ~ 2.5 x 10 5 M 0 and an age ~1 Myr, 
it is 2-4 times less massive and younger than previously estimated. It is within a 
dust cloud with Ay ~ 50 mag, and shows a clear nearIR excess, likely from hot 
dust. The second radio nebula cluster is also ~1 Myr old, confirming the extreme 
youth of the starburst region. These two clusters account for about half of the 
ionizing photon rate in the radio nebula, and will eventually supply about 2/3 of 
the mechanical energy in present-day shocks. Additional sources are required to 
supply the remaining ionizing radiation, and may include very massive stars. 

Subject headings: galaxies: general - galaxies: starburst -galaxies: dwarf - galax- 
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ies:individual (NGC5253) - galaxies: star clusters: general 


Introduction 


Local dwarf starburst galaxies are close counterparts to the high-redshift star forming 
systems that built today’s galaxies via interactions and mergers. The investigation of nearby 
dwarfs that are undergoing starburst events may, thus, shed light on the way galaxies assem¬ 
ble their stellar populations across cosmic times, and on the role young massive star clusters 
have in the energy and mechanical output of star formation. 

The extreme youth of the starburst in the center of the dwarf galaxy NGC 5253 has been 
established by many investigators (e.g., van den Bergh|l980 Moorwood fe Glass|1982 Rieke, 
Lebofsky, fe Walker||1988[ |Caldwell fe Phillips 1989 Beck et al. 1996 Calzetti et al. 1997 


Pcllcrin & Robert] 2007), and continues to be supported by recent data. The majority of the 
star clusters located within the central ~250-300 pc has ages in the range from ~10 6 yrs to a 


few 10' yrs ( 

Calzetti et al. 

1997; 

Tremonti et al. 

2001 

Harris et al. 2004 

Chandar et al. 2005 

Cresci, Vanzi & Sauvage 

2005| 

de Grijs et al. 

2013) 

. A few older clusters, up to ~10 10 yr 


in age, are located farther away from the galaxy’s center (Harbeck, Gallagher & Crnojevic 


2012 de Grijs et al. 2013). The youth of the central starburst is further supported by the 


absence of detectable non-thermal radio emission ( 

Beck et al. 

1996 

) and the presence of 

strong signatures from Wolf-Rayet stars (Campbell, Terlevich & Melnick 

1986] 

Kobulnicky 

et al. 1997 Schaerer et al. 1997) Lopez-Sanchez et al. 2007 

VIonreal-Ibero et al. 2010 


Westmoquette et al. 2013), which set a limit of <3-4 Myr to the most recent episode of 


star formation. The age range of the diffuse UV stellar population (Tremonti et al. 2001 


Chandar et al. 2005) and the recent star formation history of NGC 5253 (McQuinn et al.|2010 


Harbeck, Gallagher & Crnojevic 2012) indicate that the star formation has been elevated, 
relative to the mean Hubble time value, for the past ~5x 10 8 yr. 

The question of how to sustain continuously elevated star formation, possibly in the 
form of subsequent bursts, in NGC 5253 has been tackled by many authors. An encounter 
with the relatively nearby grand-design spiral M83 about 1 Gyr ago has been suggested 


as one of the potential initial triggers (e.g. van den Bergh 1980 Caldwell & Phillips 1989). 


M83 is located at a distance of 4.5 Mpc (Thim et al. 2003) and is 1 0 54' to the NW of 
NGC 5253; thus, M83 is separated from NGC 5253 (at a distance of 3.15 Mpc Freedman 


et al. 2001; Davidge 2007) by about 1.35 Mpc. Although the distance is significant, and 


although Karachentsev et al. (2007) place NGC 5253 in the neighboring Gen A subgroup, 


Lopez-Sanchez et al. (2012) argue that NGC 5253 is located at the boundary between the two 
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subgroups of Cen A and M 810 A past interaction with the latter galaxy could explain the 
tidal extension in HI to the SE of M 83 and the extension to the North of the HI distribution 
in NGC5253. These tails could be providing the fuel for the past and current bursts of star 


formation in NGC5253, in the form of in-falling metal-poor HI clouds (Lopez-Sanchez et 


al. 2012). The in-falling clouds convert to higher density molecular gas once they enter the 


central galaxy region and mix with the local interstellar medium (ISM; Turner, Beck & Hurt 


1997 Meier, Turner & Beck 2002 Turner et al. 2015). The potential entrance ‘channel’ for 


the gas is defined by the only prominent dust lane, which bisects the galaxy roughly along 
the minor axis and emits in CO (Walsh & Roy 1989 Meier, Turner & Beck |2002 Turner et 


al. 2015). 


Thus, the current starburst in NGC5253 is possibly the latest episode of a series of such 
feeding events, which are still ongoing. The dust-corrected UV and Ha luminosities both 
provide a consistent value of the star formation rate, SFR=0.1-0.13 M 0 yr^ 1 (Calzetti et 


al. 2004, 2015), also in agreement with the SFR derived from the total infrared emission, 
SFR(TIR)= 0.1 M 0 yr -1 (using L r //j=3.7xl0 42 erg s _1 , which we calculate from the Spitzer 


imaging data of Dale et al. 2009). Radio measurements at 0.3, 0.7, 1.3, and 2 cm of the 


free-free emission (Turner, Beck & Ho 2000 Meier, Turner & Beck 2002 Turner & Beck 


2004) yield a SFR~0.3-0.36 M 0 yr -1 which is roughly a factor of three higher than what is 
obtained from the TIR and from the dust-attenuation-corrected UV and Ha. The relatively 

100 pc and ~160 pc, respectively (Calzetti et al. 2004), 
3.5 M 0 yr -1 kpc~ 2 , confirming the star 


small Ha and UV half-light radii, 
imply a high star formation rate density, T, S fr 


bursting nature of the galaxy (Kennicutt & Evans 2012). The specific SFR of NGC5253 
is sSRF~0.6-1.4x 10~ 9 yr -1 , for a stellar mass M* ~2.2xl0 8 M 0 (Calzetti et al. 2015); the 


galaxy lies above the Main Sequence of star formation, i.e., the SFR versus stellar mass 


relation, for local galaxies (Cook et al. 2014), as expected for a starburst. 


Most of the current activity is coincident with a centrally concentrated, dusty radio 
source about 15-20 pc in extent, which we term the ‘radio nebula’. This has enough free- 
free emission to require one or more <3 Myr old star clusters with total mass M 10 6 M 0 
(Turner, Beck & Ho |2000 Turner & Beck 2004), for a 0.1-120 M 0 Kroupa stellar Initial 


Mass Function (IMF Kroupa 2001 Leitherer et al. 1999). The ratio of the stellar mass to 
gas mass in the region suggests a star formation efficiency around 60%, or about 10 times 


Where are still uncertainties in the actual distances of both M83 and NGC5253. 


Karachentsev et al. 


(2007) place M83 at a distance of 5.2 Mpc, and Sakai et al. |2004) place NGC 5253 at a distance of 3.6 Mpc, 
the latter much closer to the distance of Cen A, 3.8 Mpc. This has led Karachentsev et al. (2007) to associate 
NGC 5253 to the Cen A subgroup. In the latter case, NGC 5253 may have interacted with Cen A, instead 
of M83, in the past. 
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higher than that of Milky Way clouds (Turner et al. 2015). At least two distinct young star 
clusters are identifiable in the region, one of which is heavily attenuated by dust, and has 


been associated with the peak of emission at 1.3, and 2 cm by Alonso-Herrero et al. (2004); 
this source has angular size 0.05"x0.1" (~0.8x1.6 pc 2 ) and is associated with ~20%-30% 


of the ionizing photons in the radio nebula (Turner & Beck 2004). The other cluster is 


also affected by the dust contained in the radio nebula, but to a much smaller degree; it is 
relatively bright in the UV, and it corresponds to the peak of observed Hct emission in the 


galaxy (Calzetti et al. 1997). Alonso-Herrero et al. (2004) associate this UV-bright cluster 
with the secondary peak of emission at 1.3 cm (Turner & Beck 2004). The radio nebula is 


driving most of the ionization in the galaxy, and the past and on-going starburst has been 
stirring the surrounding ISM, both chemically and energetically. 

NGC5253 is one of the few known cases containing regions of well-detected nitrogen 
enhancement, likely due to localized pollution from Wolf-Rayet stars in the area of the radio 


nebula ( 

Walsh & Roy|l989 

Kobulnicky et al. 1997 

Schaerer et al. 

1997; 

Monreal-Ibero et al. 

2010 

Monreal-Ibero, Walsh & Vilchez 

2012 

Westmoquette et al. 

2013 

). However, no other 


chemical ‘anomalies’ have been convincingly detected. Tentative reports of He enhancement 


(Campbell, Terlevich & Mclnick 1986; Lopez-Sanchez et al. 2007), also a potential sign of 


pollution from Wolf-Rayet or other very massive stars, have been recently cast into doubt 
(Monreal-Ibero et al. 2013). The galactocentric profile of the oxygen abundance is fairly flat 


(Westmoc 


& Vilchez 


uette et al.|2013 ), with a mean value of 12+log(0/H)=8.25 (Monreal-Ibero, Walsh 


2012), or about 35% solar^j and with some scatter depending on assumptions for 


the electron temperature zone model (Westmoquette et al. 2013). This value of the oxygen 


abundance is similar to the one reported by Bresolin (2011), 12+log(0/H)=8.20±0.03. 


The ionized gas emission shows evidence of feedback from previous activity in the galaxy: 


filaments, shells, and arches characterize the Ha distribution (Marlowe et al. 1995 Martin 


1998 Calzetti et al.||1999 ), closely followed by the X-ray emission tracing the hot gas (Strick¬ 
land & Stevens 1999 Summers et al.|2004). The Ha is mostly photo-ionized but also includes 


a non-negligible fraction, up to 15% in luminosity, of shock-ionization (Calzetti et al. 1999 
2004 Hong et al. 2013); the ~kpc-size shells expand at a velocity of ~35 km s^ 1 (Marlowe 


et al. 1995) and have ages around 10-15 Myr (Martin 1998). Thus, the clusters and stars 


located in the starburst have a major impact on a number of observable characteristics of 
this galaxy, which would otherwise appear to be a rather unremarkable early-type dwarf. 

Despite ample evidence for mechanical feedback, we will assume in this work that only a 
small fraction of ionizing photons escapes the galaxy. This is true for local starburst galaxies 


2 We adopt 12+log(O/H) 0 =8.69, for the solar oxygen abundance value (Asplund et al. 


2009) 
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in general, where escaping fractions are less than 3% (e.g., Grimes et al. 2009s Leitet et ah 


2013). Recently, Zastrow et ah (2013) have suggested that these fractions may be lower limits 


due to the presence, in several starburst galaxies, of optically thin ionization cones, which 
may act as channels for the escape of ionizing photons. These will remain mostly undetected 
due to the random orientation of the cones relative to the line of sight. In NGC5253, the 
putative ionization cone is coincident with the dust lane (Zastrow et ah 2011) and with 
optically-thick CO(3-2) emission (Turner et ah 2015). Thus, while the properties of this 
feature are consistent with photoionization by escaping radiation, the ionization cone may 
be dusty and optically thick. Furthermore, the high ionization levels that mark this feature 
in NGC5253 are found in only one direction, so the solid angle of the escape zone is likely 
to be small (Zastrow et ah 2013). In what follows, we assume negligible escape of ionizing 
photons from NGC5253, although the issue remains open. 


The new high-spatial resolution UV observations presented here provide an essential 
wavelength for probing the massive star population and the impact of dust extinction in the 
radio nebula. Our goal is to quantify the properties of the star clusters in the radio nebula, in 
order to better understand their energetics and role within the NGC5253 starburst. To this 
end, we study the stellar population content of the two star clusters using SED-modcling 
techniques on UV-optical-nearIR HST photometry. The photometric stellar continuum 
bands are supplemented with measurements of the emission lines in the light of Ha, P /3, 
and Pa, also from HST imaging, which help to further constrain the ages and masses of 
the star clusters. The robustness of the SED modeling is first tested against other bright 
stellar clusters within the starburst region of NGC5253, which are less affected by dust 
attenuation than the clusters within the radio nebula, and can, thus, provide a handle on 
potential degeneracies in the results for the latter. 

This paper is organized as follows: Section 2 describes the observations and the archival 
HST data used in this investigation; Section 3 presents the cluster selection and photome¬ 
try; Section 4 presents the synthetic photometry and the fitting approach to the observed 
one; Section 5 describes the results of the SED fitting, and provides the ages, masses, and 
extinctions of the clusters, which are further discussed in Section 6. A summary and the 
conclusions are provided in Section 7. 























2. Observations and Archival Data 


2.1. New Observations 


NGC5253 was observed with the HST Wide Field Camera 3 (WFC3) in the UVIS 
channel, in the filters F275W and F336W, on 2013-08-28, as part of the HST Treasury 
program LEGUS (Legacy ExtraGalactic UV Survey, GO-13364). A description of the survey, 
the observations, and the image processing is given in Calzetti et al. (2015). 


Briefly, the WFC3/UVIS datasets were processed through the calwf3 pipeline version 
3.1.2 once all the relevant calibration hies (bias and dark frames) for the date of observation 
were available in MAST. The calibrated, hat-fielded individual exposures were corrected for 
charge transfer efficiency (CTE) losses by using a publicly available stand-alone program^ 
These corrections were small because we used the post-hash facility^] to increase the back¬ 
ground to a level near 12 e“. The processed individual dithered images were then aligned, 
cosmic-ray cleaned, sky-subtracted, and combined at the native pixel scale using the astro- 
drizzle routine^] to an accuracy of better than 0.1 pixels. The World Coordinate System 
of the WFC3 F336W image was propagated to the other image, to obtain aligned images 
across hlters, and the images in both hlters were aligned with North up and East left. The 
images are in units of e“s _1 , which are converted to physical units using the WFC3 pho¬ 
tometric zeropoints, included as keywords in the headers of the data products and posted 
at: http://www.stsci.edu/hst/wfc3/phot_zp_lbn, The basic details of the LEGUS images for 
NGC5253 used in this paper are given in Table [lj 


2.2. Archival Images 

The HST Archive contains a rich collection of images for NGC5253. For this paper, we 
retrieved images spanning from the UV to the H-band through the Hubble Legacy Archive]^] 
(HLA), both broad and narrow-band, to cover stellar continuum as well as optical and nearIR 
emission lines. When images in similar bands were available, preference was given to those 
at the higher angular resolution (e.g., ACS/HRC images were preferred over ACS/WFC 
images). Because of the extended wavelength coverage, the images used here have been 


3 Anderson, J., 2013, http://www.stsci.edu/hst/wfc3/tools/cte_tools 

4 http://www.stsci.edu/hst/wfc3/ins_performance/CTE/ANDERSON_UVIS_POSTFLASH_EFFICACY.pdf 
5 see: http://drizzlepac.stsci.edu/ 

6 http: //hla.stsci.edu 
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obtained with different HST instruments, including the ACS/SBC, ACS/HRC, WFC3/IR, 
and NICMOS. Level 2 products were retrieved for each instrument/filter combination, imply¬ 
ing that the individual post-pipeline exposures have been aligned, cosmic-ray cleaned, and 
combined using either MultiDrizzle or ASTRODRIZZLE. The retrieved images have also been 
geometrically corrected and aligned with North up and East left. All images are provided by 
the HLA in units of e - s -1 . We convert all images, except those from NICMOS, to physical 
units using the photometric zero points appropriate for each instrument/filter combination. 
The photometric zero points of the NICMOS images are referred to the default calibration 
of the instrument in DN s _1 , thus we divide the HLA images first by the NICMOS Camera 2 
gain (5.4 e _ /DN) in order to apply the published zero points. The details for each image 
product are listed in Table [TJ 

Although the archival data display a range of depths (as indicated by the large range 
of exposure times in Table [I]), all sources we study are detected with S/N>100 in the broad 
and medium band Liters. The uncertainty in the photometry is driven by crowding and 
uncertainties in the aperture corrections, rather than S/N limitations. 


2.3. Additional Processing 

Improved alignment of all the images, both new and archival, is accomplished using the 
mAh!] tasks geomap and geotran and a sample of stellar sources in the ACS/HRC images 
as reference. The HRC images are preferred over others, because they have the smallest 
native pixel, which, in this case, drives the angular resolution of the final images. Thus, 
we elect to preserve as much as possible the highest angular resolution, even if it results in 
oversampling some of the lower resolution images. For the same reason, the aligned images 
are all re-sampled to the pixel scale of the ACS/HRC, 0.025" pix” 1 (Table [TJ. 

After alignment, all images dominated by stellar continuum, i.e., all Liters except F658N, 
F129N, and F187N, are converted to physical units of erg s -1 cm -2 A -1 using the most up- 
to-date values of the PHOTFLAM keyword as posted on the relevant webpage for each 
instrument (see example for WFC3 in previous section). 

The nebular continuum and line emission from the central radio nebula signiLcantly 
contaminate the Luxes in the broad band Liters, for the clusters both within the nebula 


'IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associ¬ 
ation of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National 
Science Foundation. 
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and in the surrounding region. For instance, the presence of line emission in the F814W 
filter increases the flux measured for our individual sources between a few percent and a 
factor >3, depending on the location of the source. This effect has been noted as a problem 


for measurements of young sources by others (Johnson et al. 1999 Reines et al. 2010). 


Contamination by emission lines of broad band filters, in turn, affects the derivation of the 
line flux intensities themselves, since the broad band images are used for the subtraction of 
the stellar continuum from narrow band images. Nebular continuum will not have the same 
effect, since it is present in both broad and narrow band filters. 

We derive emission-line-free images for the most affected among our filters: F435W, 
F814W, and F110W. The F110W filter receives most of the emission line contribution from 
P/3, and we use iterative subtraction between the F110W and F128N filters to remove the line 
contamination. For the lines affecting the F435W and F814W filters we do not have direct 
imaging in the corresponding narrow-band filters. We thus use the 3,200-10,000 A spectrum 


of Storchi-Bergmann, Kinney & Challis (1995) of the central 10 ,, x20" region of NGC5253 


convolved with the F435W and F814W transmission curves to estimate the emission line 


contamination in these filters. The spectrum by Storchi-Bergmann, Kinney & Challis (1995) 


covers a sizable fraction of the region of interest here, along the E-W direction, and is thus 
representative of the excitation conditions in the center of NGC 5253. The Ha image derived 
from the F658N filter (see below) is then rescaled to the intensity of the emission lines and 
subtracted from both the F435W and F814W images. This process converges within two 
iterations. The remaining broad and medium band images are not significantly contaminated 
by emission lines, as estimated from the same spectrum. 


Emission line images are derived directly from the narrow-band filters, after subtracting 
the underlying stellar and nebular continuum. All narrow-band filters are converted to 
monochromatic fluxes (erg s -1 cm -2 A -1 ), before performing continuum subtraction. The 
continuum images are derived as follows. For the F128N image, which contains the P/3 line, 
the rescaled, nebular-line-subtracted F110W image is used. Although straightforward, this 
method can include hard-to-quantify uncertainties, if there are significant color changes in 
the stellar population across the field-of-view. For the F187N image, which contains the Pa 
line, we employ the rescaled F190N narrow-band image (Table [TJ) , which is free of emission 
lines and of any complications induced by potential color changes across the held. For both 
the FI 10W and F190N images, the rescaling factors are determined from emission-free point 
sources. For the F658N Elter, which contains Ha+[NII], we create a continuum image by 
interpolating the flux-calibrated F550M and line-emission-subtracted F814W images. The 
resulting image is then subtracted from the flux-calibrated F658N image, without rescaling. 

The line emission images are then converted to units of erg s -1 cm -2 by: (a) multi- 
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plying each image by the filter bandpas^] (72A, 159A, and 188A, for F658N, F128N, and 
F187N, respectively); and (b) correcting for the filter transmission curve values at the lo¬ 
cation of the redshifted lines. We remove the [Nil] emission from the F658N image using 
[NII](6584A)/Ha=0.084 from Moustakas & Kennicutt (2006), which is close to the value 
obtained from the spectrum of Storchi-Bergmann, Kinney & Challis (1995), and the atomic 
ratio [NII](6548A)/[NII](6584A)=0.3. 


3. Cluster Selection and Photometry 


Two cluster candidates are selected within the radio nebula (Figures la and [2]): one 
corresponding to the observed peak in Ho and the other corresponding to the observed peak 
in Pa. Measurements at 7 mm indicate a size of ~1."2 (~18 pc) for the radio nebula (Turner 


& Beck] 2004), as shown by the orange circle in the left-hand-size panel of Figure [2] The two 


peaks, Ha and Pa, are separated by about 0.46" mainly along the E-W direction (Figure [2]), 
corresponding to a spatial separation of ~7 pc, with the Pa peak emission located to the 
west of the Ha one. For each peak, the other line is also present, but not as prominently. 
The Ha peak, called ‘5’ in Figure la, has both P (3 and Pa emission spatially coincident 
with each other, and also with the continuum emission, within the accuracy that can be 
established from the image-to-image resolution differences (column 3 of Table [l]). 


The cluster candidate corresponding to the Pa peak, called ‘IP in Figure [laj is slightly 
offset, by about 0.1" (~1.5 pc), to the East of the centroid of the Ha emission closest to 
it, while the P/3 centroid falls in-between the peak locations of the other two line^} This 
gradual transition as a function of increasing wavelength suggests that the offsets between 
the peaks of the hydrogen emission lines are likely due to variations in the dust optical depth, 
rather than the presence of separate sources of emission. Although the latter scenario cannot 
be completely ruled out, we will assume in this work that the slightly spatially-shifted lines 
all originate from the same source. A visual inspection of the continuum images shows that 
the spatial shift occurs between the J (F110W) and H (F160W) images, and no shift is 
obviously present at shorter wavelengths; the centroid of the source in the NICMOS F110W 
image coincides with the centroids in the shorter wavelength images, while the centroid in 


s The filter bandpass is defined as the filter rectangular width, i.e., the equiva¬ 
lent width divided by the maximum throughput within the filter bandpass. See, e.g., 
http: / / www.stsci.edu/hst / wfc3/documents/handbooks / currentIHB / c07_ir06.html 

9 Centroids of local emission peaks can be determined with an accuracy of about l/5th-l/10th pixel, 
which, for the low-resolution WFC3/IR images, corresponds to a location accuracy of better than ~0.02". 
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the NICMOS F160W image coincides with the centroid of the Pa peak. 


Both clusters 5 and 11 are close to the peaks of free-free emission at cm wavelengths 
studied by Turner, Beck & Ho (2000) (Figure |2j). Cluster 11 is within 0.18", towards the 
S-W direction, of the peak at both 1.3 cm and 2 cm, while cluster 5 is ~0.18" to the South 
of the secondary peak at 1.3 cm. The coincidence between the sources would be increased 
if the relative astrometry between the HST and the cm-wavelength observations were off 
by about 0.2" along the N-S direction. This is consistent with the 0.1"-0.3" uncertainty 


of the absolute astrometry for HST images (e.g., Koekemoer et al. 2006). We thus believe 


the optical and radio peaks to be actually coincident, in agreement with the assumption 


of Alonso-Herrero et al. (2004); the observed offsets are likely due to small errors in the 


absolute reference frames of the two sets of data. 


An additional nine star clusters, all visually identified as local peaks of emission in the V 
(F550M) band and all brighter than my- =18.7 mag (My < —8.8 mag), are selected in order 
to perform tests on the SED fitting approach we adopt for this study. We ensure that the 
selected sources are clusters by requiring that each source’s FWHM is at least 50% broader 
than the stellar PSF. Our compilation brings the total number of star clusters investigated 
here to 11, whose locations are identified in Figure la and best-fit ages in Figure lb (see 
next section). 


Photometry is performed for all 11 clusters in multiple ways. Our default photometry 
uses an aperture of 5 pixels radius (0.125" ~1.9 pc) with the background measured in an 
annulus with inner radius of 20 pixels and 3 pixels wide. We perform visual inspection of 
the sky annuli for each cluster to ensure that they are not affected by contamination from 
surrounding bright stars/clusters. We also run tests using sky annuli with inner radius in 
the range 15-20 pixels and width in the range 3-6 pixels, in order to quantify the effects 
of background contamination. The resulting photometry varies by less than 8%, a much 
smaller uncertainty than those introduced by other effects (e.g. aperture corrections) as dis¬ 
cussed below. For the broad-band filters, we perform photometry on both the nebular-line 
subtracted and unsubtracted images, which we will compare to appropriate synthetic pho¬ 
tometry from stellar population synthesis models (with and without nebular line emission, 
see next section). For the emission lines, we perform photometry on the stellar-continuum 
subtracted images. 


As some of the clusters show a complex structure (typically elongated), we also perform 
larger-aperture photometry, with 10-15 pixels radii (and up to 20 pixels for the emission 
lines), for the more spatially-isolated clusters. We use this larger aperture photometry as 
a check for our aperture corrections, especially for the WFC3/IR images, which have pixel 
size comparable to the radius of the default photometric aperture. We choose not to adopt 
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the larger radius apertures as default for photometry, because a few of the 11 star clusters, 
including both clusters in the radio nebula, are located in crowded regions. 


The aperture corrections are determined from isolated star clusters found around the 
region where our 11 target clusters are located. We derive separate corrections for each 
instrument/filter combination. For the medium/broad-band filters (stellar continuum), 
the aperture corrections needed to bring the 5-pixels radius photometry to the infinite- 
aperture equivalent range from a factor 1.7 (WFC3/UVIS/F275W) to a factor 2.25 (NIC- 
MOS/NIC2/F160W), with a larger value, 2.81, for the WFC3/IR/F110W instrument/filter 
combination. For the emission lines, the aperture corrections for a 5-pixels radius are sig¬ 
nificantly larger, between a factor 3.7 and 7.8, which accounts for the more extended nature 
of the nebular emission. As expected, the aperture corrections decrease significantly, with 
values ranging from 5% to 20% for the stellar continuum filters, and from 40% to 70% 
for the emission lines, when a 15-pixel radius aperture is used for photometry. Compar¬ 
isons between our default aperture and larger-aperture photometry indicates uncertainties 
of ~15% for all UV-optical medium/broad-band filters, 20% for the NICMOS stellar con¬ 
tinuum filters, and 35% for the WFC3/IR/F110W filter; for the emission lines, we derive: 
~30%-55%-35% uncertainty for Ha, P (3, and Pcc, respectively. The larger uncertainty for 
the WFC3/IR photometry simply reflects the larger pixel scale of these images. The aper¬ 
ture corrections are the largest source of uncertainty for the stellar continuum filters; the 
emission lines suffer from an additional (smaller) uncertainty due to the underlying stellar 
continuum subtraction. Together with small registration offsets, this is especially a limita¬ 
tion for the P/3 photometry, despite having one of the deepest among our exposures. The 
shallow depth of the exposure is an additional limitation for the Pa image (Table [I]). The 
combination of all uncertainties, excluding the aperture correction ones, gives la depths of: 
L(Ha)=2.3x 10 35 erg s -1 , L(P/3)=3.5xl0 35 erg s _1 , and L(Pa)=2.8xl0 35 erg s _1 . 


All photometry is corrected for foreground Milky Way extinction, using the extinction 
curve values listed in Tableland the color excess E(B—V)=0.049 from Schlafly & Finkbcineri 
(2011, as retrieved from the NASA/IPAC Extragalactic Database). The values of Table [2] 


can be directly applied to photometric measurements only for small values of the color excess, 
typically E(B—V)<0.1, since color variations across the hlter bandpass are typically small; 
for larger values of the color excess, the extinction correction should be applied to the source’s 
SED before convolution with the telescope/instrument/filter response curve. 

Table [3] lists the 11 star clusters, their coordinates, and the cross-IDs with other studies 
(Calzetti et al. 1997 Harris et al. 2004 de Grijs et al. 2013), where availably} together 


10 


Our coordinates are slightly offset relative to those of 


de Grijs et al. 


(2013) by An = -0.057 s and 





















14 


with photometry, Ha equivalent widths (EW), and the color excess as inferred from the 
hydrogen emission line ratios. The listed photometry is for the measurements performed 
in the 5-pixels apertures, corrected for foreground Milky Way extinction and for aperture 
effects; in the case of continuum images, the photometry is from the original images, which 
include contribution from emission lines. Overall, the photometry of cluster 11 has larger 
uncertainties than that of the other clusters in our sample, due to its low flux densities, 
which are from a few times to over an order of magnitude fainter, depending on wavelength. 


The EWs of Ha, calculated from the ratio of the emission line flux to the stellar contin¬ 
uum flux density (interpolated from emission-line-subtracted images, see previous section), 
are given as a range: the smaller value corresponds to the ratio of line-to-continuum for 
measurements within the 5-pixels radius aperture; the larger value corresponds to the ratio 
obtained after both line and continuum have been corrected for aperture effects. The color ex¬ 
cess values are derived from the line ratios Ha/P/3 and Ha/P a using the selective extinction 
values that can be derived from Table |2j i.e., k(Ha)—k(P/3)=1.70 and k(Ha)— k(Pa)=2.08, 
and the simple assumption of a foreground dust screen. For the intrinsic line ratios we adopt 
Ha/P/3=17.57 and Ha/Pa=8.64, which are appropriate for HII regions with electron tem¬ 
perature T e rv_/ 11,500 K, measured for NGC5253 ( |Lopez-Sanchez et al.||2007j ) . We do not 
report the color excess derived from the ratio P/3/Pa, since the selective extinction between 
the wavelengths of these two lines is small, and thus the resulting colors excess is subject to 
large uncertainties. 


Two sets of values correspond to measurements performed at similar or close wave¬ 
lengths, but with different instruments (Tablc[TJ) : the WFC3/UVIS/F336W and the WFC3/IR/F110W 
measurements can be compared with the ACS/HRC/F330W and NICMOS/NIC2/F110W 
measurements, respectively. A close inspection of the photometry listed in Table [3] shows 
that the photometry in the two blue filters, WFC3/UVIS/F336W and ACS/HRC/F330W, is 
usually comparable to better than 15% (~0.07 in log scale), with the exception of cluster 11, 
where the difference is about 25% (~0.1 in log scale). We attribute the discrepancy to the 
difficulty of determining the background level around this highly obscured star cluster; how¬ 
ever, even in this case the difference in photometry is still within the combined 1 a error of 
the two measurements. Conversely, the photometric values in WFC3/IR/F110W and NIC- 
MOS/NIC2/F110W tend to be more discrepant with each other, with differences that range 
from 10% to 40% (0.04 to 0.15 in log scale). There is no obvious trend for one measurement 
to be systematically higher or lower than the other, although the NICMOS/NIC2/F110W 
measurement is more frequently the lower value. As the NICMOS /NIC2/F110W filter is 


= -0.15". 
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at slightly shorter wavelength than the WFC3/IR/F110W, its photometry values should be 
higher, thus the observed discrepancy is likely a combination of measurement uncertainties 
and, possibly, some systematic calibration offset. Similarly to the other pair of filters, the 
discrepancies are within the combined 1 a error of the two measurements. 

Color-color plots of the 11 clusters in selected bands are shown in Figure [3| together 
with the tracks of model stellar populations (section 4). These plots are only shown to guide 
intuition, and will not be used to derive the physical properties of the star clusters. 


4. Synthetic Photometry and Fitting Approach 

Spectral energy distributions (SEDs) from the UV to the nearIR are generated using 


the Starburst99 (Leitherer et al. 1999, version as available in early 2014) spectral synthesis 


models, using instantaneous star formation, with a Kroupa (2001) IMF in the range 0.1 


120 M 0 and metallicity Z=0.004 (~30% solar), which is the closest value to the measured 
oxygen abundance of NGC5253 and for which models are available. We produce models 
using both the Padova with AGB treatment^ and the Geneva tracks (Meynet et al. 1994 


Girardi et al. 2000 Vazquez & Leitherer 2005). Since the clusters under consideration tend 


to be massive, M>10 4 M 0 , we expect minimal impact from stochastic sampling of the IMF 


(Cerviho & Luridiana 2004), and use the default deterministic models. The Starburst99 


models include nebular continuum, but not nebular emission lines. These are added by 


Yggdrasil ( 

Zackrisson et al. 

2011 

), which uses Starburst99 stellar populations as an input 

for CLOUDY ( 

Ferland et al. 

2013 

). For Yggdrasil, we adopt a 50% covering factor for 


the ionized gas, meaning that only 50% of the nebular emission is spatially coincident with 
the star cluster. This attempts to reproduce the observed trend for nebular emission to be 
more extended than the stellar continuum (see previous section). Models with and without 
emission lines are generated for the age range 1 Myr-1 Gyr in steps of 1 Myr in the 1- 
15 Myr range, 10 Myr in the 20-100 Myr range, and 100 Myr in the 200-1,000 Myr range. 
Instantaneous models are assumed here to reasonably represent the population of individual 
star clusters. 

The SEDs produced by both Starburst99 and Yggdrasil are attenuated with: a star- 


burst attenuation curve (Calzetti et al. 2000), and a Milky Way, an LMC and an SMC 


extinction curve (as parametrized by Fitzpatrick 1999). For the extinction curves, we adopt 


11 Our clusters are young enough, <15 Myr, that use of the Padova tracks without AGB treatment yields 
identical results. 
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a foreground dust geometry (Calzetti 2001) of the form: 


F(\) out = F{\) mode M-^E{B-V)m]^ 


( 1 ) 


and both cases of equal and differential attenuation for the nebular gas and stellar continuum; 
for the differential attenuation, we assume that the stellar continuum is subject to half the 


attenuation of the nebular gas (Calzetti, Kinney & Storchi-Bergmann 1994; Kreckel et ah 


2013). For the case of the starburst attenuation curve, the dust geometry is ‘built-in’ into 


the functional form of the curve, and the differential attenuation between gas and stars is 
part of the way the curve itself was derived. We thus end up with seven different models 
for the dust attenuation: one attenuation curve and six extinction curves (three times two 
different ways of attenuating gas and stars). We generate the models in the color excess 
range E(B—V)=0-3 mag, with step 0.01. 

We will see in the next section that cluster 11 cannot be easily explained by fore¬ 
ground extinction/attenuation only. For this case, we generate models in which the dust 
and stars/gas are uniformly mixed together, according to the formula: 


F(X)out = F{\) model {\ - e [- a921£(B - y)fc(A)1 }/[0.921K(S - V)k{ A)], 


( 2 ) 


with the color excess E(B—V) in the range 0-20 mag. Although the uniformly mixed geome¬ 
try is likely to be an oversimplification of the complex environment surrounding cluster 11, it 
helps explain many of the properties of the star cluster. Throughout this paper, we will call 
‘front-to-back optical depth’ the quantity Ay=3.1 E(B — V) k(V) from the mixed geometry. 

The dust-attenuated SEDs are then convolved with the transmission curve of the filter 
plus the HST optics to produce synthetic luminosities, that are normalized to the default 
mass of Starburst99, 10 6 M e . 

We use ^-minimization between the models and the data, taking into account the 
measurement uncertainties, to obtain the distribution of solutions and the reduced y 2 value 
for each. We then plot the distribution of solutions within the 99% significance level for the 
appropriate number of degrees of freedom, and select the best values and the uncertainty 
for the age, color excess, and mass of each star cluster based on the shape of the reduced 
y 2 probability distribution. We fit only the stellar continuum (medium/broad-band filter) 
photometry up to and not including the H-band. Both the J-band and the ff-band can 


be heavily affected by the presence of small numbers of red supergiant stars (Cervino & 


Luridiana 2004 Gazak et al. 2013 de Grijs et al. 2013). In order to retain as much as 


possible of the wavelength baseline, we include the J-band in our fits, but exclude the Id- 
band, and we only use it as a sanity check on our results. We use the hydrogen emission line 
intensities and the Ha EW as a check on our solutions, by deriving an approximate age from 























17 


the Ha EW and a range of color excesses from the emission line ratios. We do not include 
the emission lines in the fit directly, since these can be affected by feedback effects from the 
star clusters (e.g., supernova explosions, which begin within the first 3 Myr, can eject gas 
from the cluster’s surroundings and lead to an underestimate of the emission line intensity), 
especially for the massive clusters we are studying. 

As presented in the next section, some of the star clusters have best fit ages around 
1 Myr. This implies that pre-Main-Sequence stars could be present and contribute to the 
observed SEDs. Our models do not include pre-Main-Sequence stars, and this should be 
taken as a limitation to our approach. 

We derive three solutions for the age, color excess, and mass of each star cluster from 
SED fitting. Two are based on the full wavelength coverage from rv./ 1,500 A to -11,000 A 
(7 data points = 3 degrees of freedom, we average together the two measurements in U 
and the two measurements in J, to produce one single data point at each wavelength), 
using Starburst99 and Yggdrasil models for the nebular-line-subtracted and unsubtracted 
photometry, respectively. The solutions from the comparison of the unsubtracted photometry 
with the Yggdrasil models are our reference values. We use the sets of solutions from 
the subtracted photometry plus Starburst99 models as a comparison, in order to evaluate 
how well CLOUDY reproduces the conditions of the nebular gas in each star cluster. This 
is particularly important for the central clusters in NGC5253, where the strong ionized 
gas emission can affect the measurements (e.g., by leaving residual emission in the stellar 
continuum bands). A third set of solutions is based on using only 5 bands (F275W, F336W, 
F435W, F550M, and F814W = 1 degree of freedom) for the best fits. This third set enables 
us to compare the solutions obtained from the more restricted wavelength range (which is 
the common situation for galaxies in the LEGUS and other projects) against those, possibly 
more secure, obtained from the more extended wavelength coverage. 


5. The Ages, Masses, and Extinctions of Bright Star Clusters 
5.1. Clusters outside the Radio Nebula 

Clusters 1-4 and 6-10 are located outside the radio nebula, although still within the 


starburst region. All except for cluster 4 have been investigated before by Calzetti et ah 
(1997), Tremonti et al. (2001), Harris et ah (2004), and de Grijs et ah (2013). All are 


younger than 15-20 Myr, as determined by those authors, using either lower resolution HST 
data, from the WFPC2, or UV spectroscopy, or a combination of ACS/HRC, WFPC2, and 
NICMOS data. Those earlier papers using broad and narrow-band photometry employ a 
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more restricted wavelength range, and generally only one emission line (Ha). In onr case, 
the availability of filters further in the UV (F125LP and F275W) provides better leverage for 
constraining ages of the star clusters from photometry, and the presence of multiple emission 
lines enables additional considerations on the physical conditions surrounding the clusters. 

The best fit ages, masses, and color excesses, with their la uncertainties, are listed in 
Table [3] for these cluster^} For each cluster, we generate separate hies sorted by reduced 
y 2 values, and listing ages, color excesses, and masses for different combinations of stellar 
tracks (Geneva, Padova) and extinction/attenuation curves (MW, LMC, and SMC, both 
with and without differential treatment of lines and stellar continuum, and starburst curve). 
These hies are used to determine both the best hts and the 99% confidence histograms, an 
example of which is given in Figure [4] for cluster 1. The histograms enable us to evaluate the 
uncertainties associated with each parameter, and these are the la uncertainties reported in 
Table [3J but do not carry information on the best hts (i.e., on which of the 14 combinations 
of stellar tracks and extinction curves provides the best ht to the measured photometry). 
We infer the best ht values by extracting the model with the smallest y 2 value directly from 
the hies, and the resulting synthetic SEDs and photometry are shown for all nine clusters in 
Figures [3] (top-left panel), [5j and [6j 

A few common characteristics emerge for all nine clusters from the exercise above. All 
are better ht by Padova stellar tracks, and, within the limit of validity of our foreground 
dust extinction assumptions, by the differential LMC or by the starburst attenuation curve. 
In this context, ‘better’ means that the reduced y 2 is at least 50%, and often more than a 
factor of 2, smaller than for all other solutions. For ages <6 Myr, the Padova tracks cluster 
around 5 Myr, while the Geneva tracks tend to cluster around 3 Myr for the best-ht values. 
There is also a transition for the best fitting dust extinction/attenuation: younger clusters 
(<6 Myr) prefer the differential LMC extinction, while older clusters prefer the starburst 
attenuation curve, which has the differential treatment of lines and stellar continuum ‘built 
in’. Thus, differential extinction/attenuation is always required by the best hts solutions, 
i.e., emission lines are required to be more attenuated than the stellar continuum. In this 
case, we expect the color excess derived from line ratios to be larger than that derived for the 
stellar continuum from SED fitting. To test this, Table [3] lists side-by-side E(B—V) values 
from SED htting and from emission line ratios (columns 5 and 6, respectively). The two 
sets of values are generally consistent with each other and, within the large error bars of the 
line-derived E(B—V), we cannot exclude that the latter can be larger than the SED-derived 


12 The masses of all clusters would increase by about 60% if NGC 5253 were located at a distance of 4 Mpc, 
instead of our adopted 3.15 Mpc. Changing the stellar IMF from Kroupa to Salpeter also increases masses 
by a factor 1.6, for the same 0.1-120 M 0 stellar range. 
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E(B—V). Indeed, Monreal-Ibero et al. (2010) finds evidence for differential extinction in 
NGC5253, with the stars being less attenuated than the gas by a factor 0.33. 


The main effect of differential extinction/attenuation between lines and continuum in 
the SED fits is to reduce the contribution of emission lines to the synthetic photometry in the 
broad/medium band filters, more than what is already accomplished by constructing models 
that assume only half of the ionized gas is in front of the clusters. A similar reduction effect 
can be obtained if the gas covering factor is lower than 0.5; indeed, the aperture correction 
for the Ha line is a factor over 2.5 larger than that for the underlying stellar continuum, sug¬ 
gesting a covering factor around 0.4. Furthermore, a decrease in the contribution of the metal 
lines (the major contributors to the broad band filters) can be accomplished by changing 
the ionization parameter in the CLOUDY models. Thus, differential extinction/attenuation 
should not be considered a unique solution in this case. 

For all clusters, we also show the NICMOS/NIC2/F160W photometry values predicted 
by the best-fit SEDs in Figures [4] (top-left panel), [5j and[6j In all nine cases, the prediction 
is within 2 a of the observational value, lending further support to our results. 

As a comparison, we report in Table [4] the best-fit ages and la uncertainties as ob¬ 
tained from fitting the photometry from nebular-line-subtracted images with Starburst99 
population synthesis models. We use a method similar to the one used for the Yggdrasil 
models to derive ages and uncertainties, with the only change that we do not need to ap¬ 
ply differential extinction, since the lines are no longer included in the SEDs (the nebular 
continuum is generally a much smaller contribution than the lines). As before, the younger 
clusters, 1 through 4, are better fit by an LMC extinction curve with Padova+AGB tracks, 
although the Geneva tracks give an almost as good best fit in all cases; the Padova and 
Geneva tracks yield a peak age of 5 Myr and 3 Myr, respectively, which implies differences 
in the best fit masses of roughly a factor of 2 (with the masses from the Geneva tracks being 
the smaller of the two). Older clusters (6-10) are better fit by Padova stellar tracks with the 
starburst obscuration curve in the nebular-line-subtracted case, as well. The masses and 
color excesses are also in agreement between the fits performed on the photometry with and 
without nebular emission line: they are well within the la uncertainty for E(B—V), and are 
within 70% of each other for the mass. The main exception is cluster 6, with an estimated 
mass that is a little over two times larger for the nebular-line-subtracted photometry than 
for the unsubtracted photometry. This cluster shows a more marked tail towards older ages 
than other clusters, which accounts for the discrepancy in the most likely mass value. 

The SED-derived ages are listed in Table [4] side-by-side with the ages inferred from 
the EW(Ha). The latter age is a short-hand for checking whether the ionizing photon rates 
recovered through the emission line are consistent with those expected from the best-fit SED. 
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For clusters ~10 Myr old and older, i.e., clusters 6-10, any agreement or discrepancy may 
be caused by uncertainties, since the Ha emission is expected to be low in this age range. It 
is, however, encouraging that, for the most part, there are no major discrepancies between 
the SED-derived and the EW(Ha)-derived ages for these clusters. For the younger clusters, 
1, 2, and 4, the two ages are in good agreement. This result also suggests that ionizing 
photon leakage outside the HII regions and direct dust absorption of ionizing photons are 
not significant within/around these star clusters. For cluster 3, the SED-derived age is 
~5 Myr, but the EW(Ha) suggests values in excess of 7 Myr. The disagreement between the 
two age indicators can be ascribed to the discrepancy between the observed and the best-fit 
photometry in the V and I bands (top-left panel of Figure [5]); the stellar continuum used 
to derive the EW(Ha) is obtained from the interpolation between these two bands. The 
observed continuum values are higher than the predicted ones, possibly due to untreated 
uncertainties in the photometry and in the nebular emission, and yield an underestimated 
EW(Ha). Within the uncertainties, no discrepancy is observed between the predicted and 
the measured, attenuation-corrected Ha luminosity (see section 5.2). 


Cluster 3 is the one with the smallest mass among the nine analyzed so far, and, at a 
mean value around 4,500 M 0 , may be showing some effects of stochastic sampling of the IMF 
( Cervino &; Luridiana||2004 Gazak et al. 2013 de Grijs et al.||2013 ). For instance, stochastic 
effects may account for the low observed photometry in the nearIR relative to expectations 
from the best fit SED model (Figure [5]): the observed SED is possibly deficient in red 
supergiants, which are bright, but rare, stars that become prominent nearIR contributors 
at approximately the age of cluster 3 (Popescu & Hanson 2010 Anders et al. 2013). The 
low mass of cluster 3 requires further verihcation that we are dealing with a bona-fide star 
cluster, rather than a single, isolated star. From its attenuation-corrected SED, we infer 
that this cluster contains at least 10 O stars, thus it is unlikely to be an isolated massive 
star. 


Fitting five bands (NUV, U, B, V, I) recovers physical parameters that are close to those 
obtained with the seven bands fits for 7/9 clusters. In all cases, the extinction/attenuation 
model needs to be specified a priori, in order to converge. We run tests using both the star- 
burst curve and the LMC extinction curve with differential treatment of lines and continuum 
for the clusters younger than 6 Myr and the starburst curve only for older star clusters. The 
goal is to check whether the starburst curve can be adopted in all cases. Figure [7] shows the 
results for Clusters 1 and 10, that are representative of the seven clusters with consistent 
solutions for seven and five bands fits. Cluster 1 displays a preference for young ages inde¬ 
pendently of the extinction/attenuation model selected (either LMC or starburst), although 
the LMC yields better reduced y 2 values, by about 50%. Cluster 10 shows the same double 
peak, with a range from ~7 Myr to 15 Myr in both cases. 
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The remaining two clusters, cluster 4 and cluster 6, yield results that are not as robust as 
the other seven clusters, when using five band fits. The left panels of Figure [8] show the 99% 
confidence histograms for cluster 4, for the two cases of seven bands (red) and five bands fits 
(black). When five bands are used, the young age of cluster 4 is only marginally preferred 
for the Yggdrasil fits (top-left panel of Figure [8]), i.e. when fitting broad/medium band 
photometry that include emission lines, ffowever, when using Starburst99 to fit photometric 
data from which emission lines have been subtracted, the preferred age for cluster 4 is 
markedly young, around 4-5 Myr, for both seven band and five band fits (bottom-left 
panel of Figure [8]). Results are more complicated for cluster 6: fitting five bands yields an 
extremely young age («l-2 Myr) when emission lines are included in the photometry (top- 
right panel of Figure [8]), but marginally prefers the older age of ~10 Myr when emission 
lines are subtracted (bottom-right panel of Figure [8]). The older age is preferred by the 
seven band photometry in both cases. For both cluster 4 and cluster 6, removing emission 
lines from the photometric data yields better agreement in the recovered ages between the 
restricted and the more extended wavelength ranges. This indicates that the method by 
which emission lines are included in Yggdrasil may not reflect reality in a small fraction of 
the sources. Verifying whether the observed effect is more general than indicated by our 
limited analysis will require a larger sample of star clusters. 


5.2. Clusters within the Radio Nebula 

The results from the previous section suggest that the SED fitting yields results that 
are internally consistent, and can be used to investigate the properties of clusters 5 and 11. 
These clusters are highly attenuated by dust. The hydrogen line ratios reveal significant dust 
effects, with foreground values Ay ~1.5 mag and 4.7 mag for clusters 5 and 11, respectively. 

The results of SED fitting for these two clusters are listed in Table [5j and shown in 
Figure [9] Both clusters are extremely young, with a best fit age of 1 Myr, which agrees with 
the age inferred from the EW(Ha). Both are better represented by Geneva stellar tracks, 
and by the starburst attenuation curve for the foreground dust. As for the other clusters, 
we have verified that the best fit solution also reproduces the intensity and ratio of the line 
emission at Ha, P/3, and Pa. This is shown in Table [6j where we list, for all clusters, the 
intrinsic Ha luminosity, obtained from the attenuation-corrected measurements of Table [3j 
and the SED-predicted Ha luminosities, obtained from the models and the best-fit ages and 
masses of Tables [4] and [5j Within the 1 a uncertainties, there is general agreement between 
the two sets of values. Major disagreements are discussed below (cluster 5) and in section 6.1 
(cluster 9). 
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For cluster 5 the use of a simple dust geometry, i.e., foreground dust, is sufficient to 
recover an excellent fit for the SED, with a reduced y 2 ~1. The color excess derived from 
the SED fitting agrees with the one derived from both line ratios: Ha/P/3 and Ha/P a. 
The emission lines are strong enough that, although the formal uncertainties are significant 
(Table [5]), both ratios yield similar extinction values. Despite the goodness of fit, not all 
observational data points match the model SED: a significant deviation (~2cr) is evident 
for the V-band (F550M), which may be due to difficulties in measuring the background 
surrounding the cluster. The observed photometry is below the model’s value, which explains 
the unphysically large EW(Ha) (Table [ 3 ]). 


The mass of cluster 5 derived from the SED fit is ~7.5xl0 4 M 0 , confirming that this 
cluster is massive enough to drive the observed ionization. The Ha luminosity predicted 
for cluster 5 is L(Ha)~5.7xl0 39 erg s -1 , to be compared with the attenuation-corrected 


from the free-free measurement of 
free-free measure likely provides 

Turner, Beck & Ho 
an underestimate 0 : 

(2000), L(Ha)~2-3x 10 39 ergs 1 (the 
the ionizing photon flux, since it is 

affected by self-absorption, see 

Meier, Turner & Beck 

2002 

). The discrepancy between the 


SED-predicted and the observed, attenuation-corrected Ha luminosity can be interpreted 
as due to either leakage of ionizing photons outside of the Hll region or direct absorption 
of ionizing photons by dust. Leakage is supported by the presence of extended ionization 
around cluster 5, which suggests that ionizing photons from this cluster reach further out 
than what we recover within our photometric apertures. We infer that 25%-50% of the 
ionizing photons produced by this star cluster leak out of the region. Direct absorption of 
ionizing photons by dust is also a potential mechanism in the dense environment surrounding 
cluster 5. However, we provide arguments both below and in section 6.6 against a significant 
contribution from dust absorption of ionizing photons in this galaxy. 


Cluster 11 requires a more complex approach to dust attenuation, in order to approxi¬ 
mate the observed SED. A simple foreground dust model is unable to provide a reasonable 
fit to the observed SED, on account of the SED being mostly flat in L(A)-versus-A. A 
combination of a homogeneous dust-star mixture and a foreground dust screen provide a 
better, albeit not perfect, fit (Figure [9]). Flat or slightly blue SEDs in the presence of sig¬ 
nificant amounts of dust tend to require models in which the dust is mixed with the stars. 
Foreground or shell models, either homogeneous or clumpy, will generally tend to produce 
too red SEDs relative to the data at blue wavelengths and too blue SEDs at red wavelengths 
(Calzetti 2001). 


The major deviations between the data of cluster 11 and the model are in the V band 
and in the J and H bands. Like in cluster 5, the V-band observations are about 2a below the 
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model, possibly an effect of background placement for the measurements. As in the other 
cluster, this accounts for the unphysically large EW(Ha) measured for cluster 11 (Table [3]). 
The other major deviation, in the nearIR bands, is in the opposite direction: the data are 
more luminous, by slightly more than 2cr, than the model expectation. A discussion of this 
deviation is deferred to section 6.3. 

Taken at face value, cluster 11 is the most massive among those analyzed here, with a 
mass M~2.5xl0 5 M 0 . The component of dust that is mixed with the stars in the cluster has 
a front-to-back optical depth of Ay ~49 mag and follows the Milky Way extinction curve; 
this curve yields a reduced y 2 that is at least two times better than any other of the extinction 
curves tested in this paper. The foreground dust component obeys the starburst attenuation 
curve, with Ay ~1.9 mag. Cluster 11 is therefore associated with about 50 magnitudes 
of optical depth in dust, although the mixed geometry allows for some of the UV-optical 
radiation to shine through. This optical depth is larger than what was determined by other 


authors. Alonso-Herrero et al. (2004) derive Ay ~17 mag for this cluster, but they only 


assume a foreground dust screen, which provides a lower limit to the actual dust optical 


depth of a region. Martin-Hernandez, Schaerer & Sauvage (2005) derive a value of Ay ~14- 


17 mag, from the 9.7 /mi silicate absorption feature; however, the relatively large slit they 
use for their spectroscopic observations (1".2 ~ 18 pc) is likely to have sampled regions of 
different optical depth, thus weighting the final result towards less extreme values of Ay. 
Calzetti et al. (1997) derive a larger range, Ay ~9-35 mag, the upper bound of which is 


more consistent with the value derived here. 

From the best-fit model SED, the predicted Ha luminosity of cluster 11 is L(Ha)~2x 10 40 erg s _1 , 
which is comparable to the measured, attenuation-corrected L(Ha) ~1.8±££ xlO 40 erg s" 1 
(Table [6j). As a reminder, the full dust model of combined mixed and foreground dust is 
applied to all three emission lines, Ha, P (3, and Pa, to derive internally-consistent values for 
their luminosities. The free-free emission within 4 pc of the radio peak corresponds to an 
intrinsic Ha luminosity L(Ha)~1.6-2.0x 10 40 erg s~ 4 (from the 7 mm measurement of Turner 


& Beck 2004). The agreement between the measured, attenuation-corrected Ha luminosity 
and the SED-predicted luminosity provides independent confirmation of the accuracy of our 
results for cluster 11. Within the uncertainty of the measurements, the mass and age of 
cluster 11 are accurate, as are the dust content and geometry; there are no major parts of 
cluster 11 that are so deeply buried in dust to be unrecoverable by our approach. It also 
suggests that, within the uncertainties, direct absorption of ionizing photons by dust does 
not appear to be a dominant mechanism at work in the area where cluster 11 resides, despite 
the large dust optical depth. 


Because of the red SED of cluster 11, a secondary best-fit solution is provided by a 
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~100 Myr old star cluster, with a mass M~3.5xl0 5 M 0 , also mixed with dust with optical 
depth Ay ~42 mag. We reject this solution on account of the strong nebular line emission in 
coincidence of the location of cluster 11. However, the uncertainties in both line emission and 
stellar continuum measurements allow the co-existence of two populations at the location 
of cluster 11: a ~2xl0 5 M e , 1 Myr old population together with a ~lxl0 5 M 0 , 100 Myr 
old population. This solution would push against the tolerance of our error bars, but is not 
formally excluded by the data. 

Other sources of uncertainty for the best-fit solutions of cluster 11 include the require¬ 
ment that emission line intensities and ratios be consistent with the model for the stellar 
continuum. If we allow the stellar continuum to be modeled independently of the emission 
lines, we obtain a larger range of degeneracies as a consequence of the larger number of 
degrees of freedom. For instance, in the absence of constraints from the emission lines, the 
dust-star mixed model produces a linear correlation between mass and front-to-back optical 
depth for Ay >15 mag: by doubling the total attenuation, one can double the cluster mass 
recovered with virtually unchanged y 2 values. 


5.3. Limitations in the Population Synthesis Models 


There are several properties of massive stars that are not currently included in the 
population synthesis models used in the present analysis. These are: stellar rotation, binary 
stars, and Very Massive Stars (VMSs). Most of these properties are still under investigation, 
and their inclusion in models is either in the early stages (e.g., rotation and binaries) or 
non existent. We provide here a brief summary of the effects we expect each to have on our 
results. 


A stellar population containing rotating massive stars will appear younger than a non¬ 
rotating counterpart, with an increase in the ionizing flux that can be as much as a factor 


two-five higher (Leitherer et al. 2014). Applied to the star clusters in our sample, such 


models are expected to yield older best-fit ages than what we derive or, alternatively, larger 
ionizing fluxes at fixed age. 

Young massive stars are predominantly found in close binaries, and can interact (e.g., 


Sana et ah 2012). The products of such interactions include massive blue stragglers formed 


through mass transfer and mergers, which may be common (de Mink 2014), and appear at 


the upper tail end of the stellar IMF (Schneider et ah 2014, 2015). The effects on young star 
clusters include apparent rejuvenation and age spreads (Eldridge & Stanway:2009). Our data 
do not include enough information to infer whether either of those effects may be present, 
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but we cannot exclude them either. 


The birth mass of stars is limited to ~120-150 M 0 in virtually all population synthesis 
models, but there is mounting evidence that stars as massive as 300 M 0 are present in 
nearby young star clusters (Crowther et ah 2010). There is still debate on whether these 
VMSs are the result of birth conditions or of mergers, and stellar tracks are being produced 
in order to further investigate this issue ( Yusof et al.|2013 Koehler et al.|2015 ). The youngest 
among our clusters, clusters 5 and 11, are also sufficiently massive that they may contain 
VMSs. Presence of VMSs increases the ionizing photon flux from a cluster. If the VMSs 
contributions to clusters 5 and 11 were similar to those found in the LMC cluster R136, the 


ionizing photon flux could be 50%-100% larger than the one currently predicted (Crowther 


et al. 

2010; 

Doran et al. 

2013) 


One final limitation of our model fitting is that the youngest synthetic population we 
consider is 1 Myr old, in line with our photometric uncertainties which yield a best accuracy 
of ~1 Myr at those young ages. The formal solutions for clusters 5 and 11 do not exclude 
that they could be younger than 1 Myr. 


6. Discussion 


6.1. Age Comparisons with the Literature 


The physical parameters of age and mass have been derived before by several authors 
for most of the star clusters in this study ( Calzetti et al. [1997 Tremonti et al. [2001 Harris et 
al. 2004; de Grijs et al. 2013). We compare our derived ages with those previous derivations 
in Table [7} We do not compare masses, as these are somewhat degenerate with ages and 
depend on additional assumptions such as the stellar IMF and the adopted distance for 
NGC5253. There is general agreement between the ages derived by all authors, with a few 
exceptions. 

The most notable discrepancies are present for cluster 9, with age estimates that range 


from 3 Myr (Tremonti et al. 2001) to 30-50 Myr (Calzetti et al. 1997). The old age derived 


by Calzetti et al. (1997) is an effect of the limited number of broad bands available to those 


authors, since they only used NUV (centered at 2250 A), V, and I ; in particular, the 
absence of a filter close to the age-sensitive U-band limits the ability to attribute ages to 
star clusters (e.g., Lee, Chandar fe Whitmore|2005 ). The younger age derived by Tremonti et 
al. (2001) is far more puzzling. These authors use UV spectroscopy for deriving cluster ages, 


leveraging the information from the photospheric lines, while all other authors derive their 
age estimates from colors or SED fitting. The UV spectrum of cluster 9 shows the presence 








































26 


of P-Cygni profiles for the lines of [NV](A1240 A) and [CIV](A1550 A), a clear sign for the 
presence of early O-stars and ages <5 Myr. However, the long slit of the HST/ST1S only 
‘grazes’ the outskirts of cluster 9, with only l/30th of the light from this cluster captured by 


the UV spectrum (Tremonti et ah 2001). Possibly, the UV spectrum is targeting a smaller 


cluster in the periphery of cluster 9. The presence of one or more interlopers is supported by 
the anomalously large (for its age and location) color excess from emission lines for cluster 9 
(Table|3]); these large values for E(B—V) are usually found in correspondence of much younger 
star clusters. Despite the large attenuation value, the emission lines in correspondence of the 
cluster are still weak for its mass (Table [6]), and there is no evidence for additional ionized 
gas in the form of shells or arcs in the region; this further supports the older age, ~10 Myr, 
for cluster 9. 


For cluster 10, the main discrepancy is given by the age reported by Calzetti et al. (1997), 


50-60 Myr, while our determination and those from other authors suggest an age around 
8-15 Myr. As in cluster 9, the discrepancy is due to the limited number of photometric 
bands available to Calzetti et al. (1997). An early suggestion that clusters 9 and 10 could 


be younger than the 30-60 Myr age inferred by Calzetti et al. (1997) came from Strickland 


& Stevens (1999), based on the measured X-ray sizes and luminosities of the super bubbles 


in NGC 5253. 

The young ages of the dusty clusters 5 and 11 have been known for quite some time 


e.g, Meurer et al. 1995 Schaerer et al. 1997 Calzetti et al. 1997), but the realization that 


the region hosts two separate clusters instead of a single one is more recent (e.g. Alonso 


Herrero et al. 2004). Alonso-Herrero et al. (2004), specifically, favors the presence of two 


star clusters over other interpretations, such as that cluster 5 may be a reflection nebula 
generated by cluster 11 (Turner & Beckj2004). Our analysis also favors the interpretation of 
two separate star clusters. The morphology of cluster 5 is consistent with that of a compact 
star cluster: a slightly-resolved (~0.1" corresponding to a deconvolved size of 1.2 pc), 
symmetric, and centrally concentrated source, similar to other clusters in the area, including 
those in the same ~1" region as cluster 5. At least an additional 6-7 fainter star cluster 
candidates are visible in the ~0.6" area surrounding cluster 5, suggesting that this one is the 
most prominent one in an association of stars and/or star clusters. 


Both clusters have been characterized as having ages ~3-3.4 Myr by Alonso-H errero 
et al. (2004), based on HST nearIR imaging and ground-based spectroscopy; the inferred 
6.6xl0 4 M 0 for cluster 5 and ~3.9xl0 5 -1.3xl0 6 M 0 for cluster 11, when 


masses are 


reported to our adopted distance and stellar IMF. For cluster 5, the difference between our 


mass and the mass derived by Alonso-Herrero et al. (2004) can be entirely attributed to 


differences in the derived ages and extinctions. For cluster 11, a major component of the 
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discrepancy is the fraction of hot dust contributing to the K-band: |Alonso-Herrero et al. 


(2004) adopt a fraction ranging from 0% to 70%, but we derive a fraction closer to 90% from 


our best fit (see section 6.3). The differences in age and dust column density and geometry 
account for the remaining portion of the discrepancy. Even with these discrepancies, the 


masses we derive for both clusters and those of Alonso-Herrero et al. (2004) agree to better 


than 60%, when using their more conservative estimate for cluster 11. 

The area surrounding clusters 5 and 11 hosts a half-dozen Wolf-Rayet stars, mainly of 
the younger WN type. This would suggest presence of stars/clusters in the narrow age range 
~2.5-3.5 Myr (Schaerer et al. 1997 Monreal-Ibero et al. 2010 |Westmoquette et al. 2013). 
This is not necessarily in contradiction with the ~1 Myr age we derive for clusters 5 and 
11. The relatively low spatial resolution of ground-based spectroscopy, ~l".5-2" (20-30 pc, 
or about the size of the panels in Figure [2J see a summary description in Westmoqu ette 


et al. 2013), has not enabled accurate location of the stars or clusters responsible for the 


Wolf-Rayet emission features. This suggests two possible scenarios. In one scenario, the WN 
stars co-exist in the region with clusters 5 and 11, and these clusters are the most recent 
‘products’ of ongoing star formation over the past few Myr. In the second scenario, the WN 


features may originate from VMSs (Crowther et al. 2010), hosted in the extremely young 
clusters 5 and 11. 


Whitmore et al. (2011) suggest a method based on the morphology of the Ha emission 


in order to classify star clusters according to their ages: the gas morphology is compact 
and coincident with the star cluster for ages <a few Myr; has a small ring-like structure in 
clusters up to ~5 Myr of age; has a large, well-formed ring surrounding the cluster for ages 


in the range 5-10 Myr; and is virtually absent in clusters older than 10 Myr. Whitmore et 


al. (2014a) adds earlier stages to this classification by including proto-cluster phases based 


on CO appearance. When evaluated according to the criteria of these papers, our clusters 
form a well defined sequence, with a close agreement between our SED-fitting ages and the 
morphological ages, and with clearly identifiable stages from 2 for cluster 11 (embedded 
cluster) to 5 for clusters 6-to-10 (intermediate/old cluster). 


6.2. Global Properties of the Clusters 


The total Ha luminosity of NGC5253, corrected for foreground Milky Way extinction 
and [Nil] contribution, but uncorrected for internal dust attenuation, is L(Ha) tofa i =2.2 x 10 40 erg s~ 


(Kennicutt et al. 2008). About 15% of the Ha is associated with shock-ionization (Hong 
et al. 2013), implying that the photo ionized Ha luminosity is L(Ha) p / lo t=1.9xl0 40 erg s 


,-i 


The sum of the observed Ha luminosity from the star clusters analyzed in this work, also 






































uncorrected for internal dust attenuation, is L(Ha) c ; usters =1.43xl0 39 erg s' 1 (from Table [ 3 ]) , 
or about 7.5% of the total Ha. Thus, the brightest 10 clusters contribute almost 10% of the 
total (observed) Ha luminosity in this galaxy; indeed, although there are almost 150 young 
star clusters in the central region of NGC5253, the vast majority tend to have low mass, 


i.e., <10 4 M 0 (de Grijs et al. 2013). 


A similar argument can be made for the FUV luminosity. The total luminosity density 
contained within the ACS/SBC frame at -—'1,500 A, corrected for foreground Milky Way 
extinction, but uncorrected for internal dust attenuation, is L(1500A)~5.8xl0 38 ergs' 1 A -1 . 
The GALEX FUV luminosity density of NGC5253 is L(1500A) tote ; ~6.6x 10 38 erg s 1 A 1 , 
i.e., only about 15% larger than the amount contained in the HST image. We assume 
this 15% discrepancy to be the upper limit to the amount of FUV light outside of the 
ACS/SBC frame, since both the SBC/F125LP filter and the GALEX/FUV filter have pivot 
wavelengths at ~1,500A. The total contribution from the 11 clusters investigated here is 
L(1500A) c / usters ~2.7x 10 37 erg s' 1 A' 1 , or about 4%-5% of the total. Star clusters contribute 
about 20% of the total observed UV light in nearby star-forming and starburst galaxies 


(Meurer et al. 1995 Maoz et al. 1996). Thus, the 11 bright clusters represents about 20%- 


25% of this contribution. 

The patchy extinction in the center of NGC5253 makes it difficult to convert the above 
numbers to intrinsic luminosities. We adopt a hybrid approach, by combining the ionizing 


photon flux measured from the free-free emission within the radio nebula (Meier, Turner & 


Beck 

2002 

Turner & Beck 

2004 

) with the virtually extinction-free Ha emission outside the 

radio nebula (Calzetti et al. 

2004 

). The two combined yield a total, attenuation-corrected 


.-1 


erg s i for the photo-ionized component, 


Ha luminosity L(H a) phot , corrected ~7.1-8.2xl0 40 
about 20%-23% of which is from outside the radio nebula. The attenuation-corrected Ha 
luminosity is, thus, ~3.5-4 times larger than the observed one for NGC5253. The SED- 
predicted Ha luminosity for all the 11 clusters is L(Ha) clusters,corrected, ~2.6xl0 40 erg s' 1 
(Table [6]), to which clusters 5 and 11 contribute 22% and 76%, respectively, and the remain¬ 
ing 9 clusters only contribute a total of 2%. These 9 clusters provide ~3% of the ionizing 
photon flux outside the radio nebula; the remaining flux is provided by the smaller clusters 


and UV-bright diffuse stars that populate the region (Hong et al. 2013). Within the radio 


nebula, clusters 5 and 11 supply about 40%-50% of the detected ionizing photon flux, i.e., 
they come short of providing the full ionizing flux by a factor ~2. We further discuss this 
discrepancy in section 6.4. 

If unimpeded by dust, cluster 11 would have an absolute magnitude My(Vega)~ —12.8, 
which is about 2-3cr above the mean of the Mf, rig toesi~SFR relation for star clusters, for 


the SFR value of NGC5253 (Whitmore et al. 2014b); this is in the bright envelope, but not 
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outside the range of statistically possible values. Indeed, the dust-free absolute V magnitude 
of cluster 11 is consistent with the results of Billett, Hunter & Elmegreen (2002), who find 
that dwarf galaxies tend to host massive clusters more frequently than do their more massive 
spiral counterparts. Billett, Hunter & Elmegreen (2002) normalize all star clusters to a 


fiducial age of 10 Myr; if cluster 11 is ‘aged’ to 10 Myr, its dust-free absolute magnitude 
would be My ~ —12.6. 


6.3. The Near Infrared Excess of Cluster 11 


The heavily dust-attenuated cluster 11 is anomalously bright in the J and H bands, 
about a factor 2-2.5 than what predicted from the best fit stellar population models. There 
are four possible causes for this, which we will analyze in turn: (1) excess nebular emission 
relative to our default model assumption; (2) red supergiants; (3) Pre-Main-Sequence stars; 
and (4) hot dust emission. Option 4 has been already considered by other authors for this 
region, and previously seen in other embedded clusters (Johnson et al. 2004), but will be 
re-analyzed here in light of our model for the stellar and dust mixture of this cluster. 


For option 1, we re-run our best-fit programs using models that implement nebular 
continuum and line emission with 100% covering factor, to mimic a tightly confined HII 
region around a star cluster. In all cases, we find that the best fit solution has a reduced 
X 2 that is at least twice as worse as the cases with 50% covering factor. The reason is 
because the nebular continuum and line emission contribute to the optical bands as well as 
the infrared ones, thus requiring a higher degree of dust attenuation to produce the observed 
SED shape. The higher attenuation pushes the UV model further away from the data, while 
not compensating enough for the high value of the nearIR data, and ultimately resulting in a 
poorer Et. Thus, a 50% covering factor is in better agreement with the data even in this more 
extreme case, and is in agreement with our measurement procedure for the stellar continuum 
and emission lines (in which we have applied larger values of the aperture corrections to the 
lines than to the continuum). 

The presence of red supergiants (option 2) is attractive because these contribute signif¬ 
icantly in the J and H bands. However, cluster 11 is too young to include red supergiants, 
which would need to be originating from another cluster/location, along the line-of-sight 
of cluster 11. These red supergiants would also need to be associated with a star cluster 
that is otherwise heavily obscured at wavelengths shorter than J and does not emit ionizing 
photons, since the entire P/3 and Pa luminosities are fully accounted for by cluster 11. The 
100 Myr old population in the two-populations solution discussed in section 5.2 is too old to 
still contain red supergiant stars. In summary, we disfavor this scenario, although we cannot 
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exclude it completely. 


Pre-Main-Sequence stars (option 3) are likely to be present in a young cluster, and to 
provide a significant contribution to the infrared emission. However, the effect corresponds 
to an increase of ~30% in the K-to-V luminosity ratio relative to a population without 
Pre-Main-Sequence stars (Zackrisson et al. 2001); this is much smaller than that the order- 
of-magnitude increase observed (see below). 

Option 4, the presence of hot dust surrounding the Pa peak, has been studied already 


by previous authors, including Vanzi & Sauvage (2004) and Alonso-Herrero et al.| (2004). 


Vanzi & Sauvage (2004) used Adaptive-Optics near infrared observations at Ks (2.16 /mi) 
and L’ (3.78 //rn) to determine and measure a peak of emission in these two bands in corre¬ 
spondence of the location of Clusters 5 and 11. The observations have FWHM=0.4 7/ , which 


is comparable to the separation between the two clusters. Vanzi & Sauvage (2004) modeled 


the SED of cluster 5 at optical wavelengths combined with their and other measurements 
at longer wavelengths, employing a physically-based model for the dust emission. Vanzi & 


Sauvage (2004), however, did not have measurements at wavelengths between I and Ks, and 


were not aware of the presence of Cluster 11, so the accuracy of their results is difficult to 
assess. Since cluster 5 does not show any evidence for near infrared excess based on our SED 
modeling, we assign those authors’ Ks and L’ emission to cluster 11. 


Alonso-Herrero et al. (2004) measured the emission from the entire central region, tar¬ 


geting an area of about 23-24 pc (l.b"xl.Q", larger at wavelengths longer than 5 /jrn), 
including both Clusters 5 and 11. They also used a physically-based model to account for 
both the stellar and dust emission of this region, but their model had difficulties in account¬ 
ing for the emission below the H band. They ascribed this difficulty to the complexity of 
the stellar and dust geometry in the region. 

Since we are interested mainly in explaining the excess emission at J and H for cluster 11, 


we use the measurements in the Ks and L’ bands from Vanzi & Sauvage (2004, their Table 3), 


which are obtained with a smaller aperture size than those of Alonso-Herrero et al. (2004). 


We attempt a simple modeling of the emission of the dust surrounding this cluster. For 
the same reason, we do not include measurements at longer wavelengths, because the lower 
resolution requires larger apertures, and more contamination from the region surrounding 
cluster 11. Using our distance for NGC5253, we derive L(2.16 /im)=2.0xl0 35 erg s~ 3 A -1 
and L(3.78 / um)=1.9x 10 36 erg s~ 3 A -1 . An extrapolation of our best-fit stellar population 
(plus dust attenuation) model shows that the stellar contribution to the 2.16 //m emission is 
about 10%, and is insignificant in the longer wavelength band. 


Adding a two-component dust emission model to the stcllar+dust attenuation model 
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reproduces the data from 1,500 A to ~4 fim reasonably well, as shown in Figure dO 


The two-component dust model consists of two modified black-bodies with temperatures 
i ~1,100 K and T^ 2 ~440 k, and emissivity e=1.8. The dust emission model parameters 
are not the result of a rigorous fitting procedure, thus the derived temperatures are approx¬ 
imate. Furthermore, an emissivity with power law of 1.8 is a gross approximation below 
20 /ini (Draine 2003), but changing the dust emissivity has a small effect on the derived 
temperatures, for our restricted wavelength range. Our two values for the dust temperature, 
1,100 K and 440 K, bracket the value of T=570 K derived by Vanzi & Sauvage (2004) for 
the inner shell of their cocoon dust model. Additionally, we have measurements in both the 
J and H bands, which further enable us to place constraints on the dust emission in this 
range and derive a higher temperature of 1,100 K. 


There is still some unaccounted for excess in the J-band, although we cannot assess its 
magnitude, given the size of our uncertainties. Thus, emission from hot dust surrounding 
cluster 11 appears to be capable of explaining most of the observed excess in the J-band 
and beyond for the SED of this star cluster. A more physically-based model for the dust 
emission would be required to infer additional properties for the dust, but this is beyond the 
scope of the present work. 


We infer the dust mass associated with these two temperature components, as a sanity 
check to the reasonability of our dust SED fit. This mass will be a lower limit to the actual 
dust mass associated with those components, since we are in the optically thick regime: the 
dust optical depth along the direction of cluster 11 is Td (4 fi m)>1.2 , which we derive using 
the gas column density discussed in the next section and the dust emissivity value at ~4 /mi 
from Draine (2004); we use the LMC dust emissivity curve from this paper, on account of 
the low metallicity in NGC5253. We derive M d us t{4 /xm)>10 M 0 , more than 99% of which 
is contributed by the cooler of the two dust components, T^ 2 ~440 K. This mass is a tiny 
fraction of the total dust mass in this region, ~10 5 M 0 , as derived by probing the cooler 
dust SED with a peak around 30-40 /irn (Vanzi & Sauvage 2004). It is, however, consistent 
with the the expectation that only the dust located in close proximity of the massive stars 
will be heated to high temperatures. For reference, the amount of hot dust surrounding 
Ultracompact HII regions in the Milky Way is calculated to be around a few to a few tens of 
M 0 (Walsh et al. 1999), consistent with our lower limit for the mass in hot dust associated 
with cluster 11. 


The presence of dust as hot as T d «1,100 K in proximity of cluster 11 supports our 
SED fitting approach of using a mixed geometry for the stars and dust. In this type of 
geometry, the dust is likely located close to the UV-emitting stars needed to heat it to 
high temperatures. Albeit large, the value of the dust temperature is still below the dust 
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sublimation temperature of ^2,000 K (Kobayashi et al. 2011), and is not dissimilar from 
that found for Ultracompact HI1 regions in the Milky Way, where single massive stars can 
heat a shell of dust in the natal cloud up to ~1,500-1,800 K (Walsh et ah 1999). 


6.4. Energy Balance Within the Radio Nebula 


Clusters 5 and 11 are the most prominent ones in the the radio nebula. These clusters 


have strong hydrogen emission line intensity and concentrated free-free emission (Alonso- 


Herrero et al. 2004 Turner & Beck 2004), although several other fainter cluster or massive 


star candidates are visible in the immediate area (<0.6" ~9 pc from cluster 5). 


The two clusters share a very young age, ~1 Myr, and a common shell of foreground 
dust, with optical depth Ay ~1.9 mag. However, they also show significant differences. 
Cluster 5 is about 3.5 times less massive than cluster 11, and is not as deeply buried in 
dust. Its location in a less dust-enshrouded environment than cluster 11 accounts for its 
prominence as the Ha peak emitter in the galaxy, and for the possibility that 25%-50% of 
its ionizing photons may be leaking out of the surrounding HII region (section 5.2). 

Cluster 11, with a best-fit mass of 2.5xlO 5 M 0 , is the behemoth in this dwarf galaxy, 
although not as extreme as inferred in earlier estimates, where values as large as rv./ 10 6 M 0 
have been suggested. While more massive than cluster 5, cluster 11 is significantly fainter 
because it is immersed in a dust cloud with front-to-back optical depth Ay ~49 mag. 


Clusters 5 and 11 contribute about 50% of the ionization in the radio nebula and about 
35% of the total in the galaxy, with most of it coming from cluster 11. While they provide 
a significant fraction of the ionizing photons in both the radio nebula and in NGC5253, 
they do not provide the totality. This poses a potential issue of energy balance, which we 
discuss in the context of the radio nebula. Within this ~20 pc region (Figure |2j) , clusters 5 
and 11 account for 50% of the free-free emission, for almost the totality of the observed Ha 
emission, but only for 23% of the observed Pa emission. The rest of the Pa emission is 
spread throughout the region of the radio nebula and has a diffuse morphology, as already 
observed by Alonso-Herrero et al. (2004). 


If we double the mass of cluster 11 in order to compensate for the factor ~2 discrepancy 
in the free-free emission, the observed stellar continuum photometry is reproduced with a 
front-to-back optical depth A y ~100 mag for the dust cloud in which the star cluster is 
immersed. The Ha emission requires the same optical depth, but the Pa emission only 
needs Ay ~20 mag, in order to recover the measured luminosity. The discrepancy in the 
dust optical depth of the two emission lines suggests that this is not a viable scenario. 
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Self-consistent solutions for the available data using a ~5xl0 5 M 0 star cluster can only 
be obtained by modeling independently the stellar continuum and the ionized emission, as 
discussed in section 5.2. We disfavor such solutions, on the basis that the lines and continuum 
of cluster 11 appear spatially correlated. 


Our original solution, with a cluster mass of 2.5 x 10 5 M 0 and front-to-back Ay ~50 mag, 
can close the gap between predicted and measured ionizing photon flux in the radio nebula, 
if young stellar populations produce more ionizing photons than what accounted for by the 


models we use. Models that include VMSs (Crowther et al. 2010 Doran et al. 2013) have 


the potential to increase the ionizing photon flux by up to a factor of 2, with smaller impact 
on the stellar continuum luminosities. The amount of diffuse Pa in the radio nebula suggests 
that about half of the ionizing photons leak out of the clusters into the region. This is similar 
to the fraction of ionizing photons leakage found by Johnson, Hunt & Reines (2009) for the 
compact, dusty clusters in the starburst galaxy SBS 0335-052. The observed Ha/Pa/free- 
free intensity ratios outside the two clusters can be fully explained by foreground dust with 
color excess E(B—V)~2 mag, a factor almost three lower than the total extinction in clus¬ 
ter 11. The gas outflow present in the area, as inferred from the broad component in the 
Ha emission (Monreal-Ibero et al. 2010 Westmoquette et al. 2013), may create favorable 
conditions for leakage of ionizing photons, by causing the ISM to become porous. 


The diffuse nature of the Pa emission outside of clusters 5 and 11 further suggests that 
this emission may originate from other stars and/or clusters in the region, some of which 
may be themselves buried in dust and undetected except at wavelengths longer than a few 
micron. 


6.5. The Environment of the Radio Nebula 


As already indicated by a number of previous studies (see Introduction), the radio 
nebula is the youngest and most active area of star formation in NGC5253. It is also the 
region of intersection between the central starburst and the dust lane, which accounts for its 
significant dust content. 


The dust cloud enshrouding cluster 11 may appear exceptional, with its Ay=49 mag, 
but it is comparable to clouds in other galaxies, including those in the Milky Way (see, e.g., 
the Ultracompact HII region W3(OH), Turner & Welch 1984). Adopting the total hydrogen 
column density-to-color excess relation of Boldin, Savage & Drake (1978), rescaled to the 
lower metallicity of NGC5253, the optical depth of the cloud corresponds to a hydrogen col¬ 
umn density N(H)~2xl0 23 cm -2 . Although large, this value is comparable to that observed 
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in some of the massive dense clouds towards the center of the Milky Way (e.g. Kauffmann 


et al. 2013). If distributed uniformly throughout a region of about 1 pc size, the observed 

g cm" 


gas column density would correspond to a dust density pd ~7xl0 22 g cm 3 , only a factor 


of a few lower than the dust densities typical of Ultracompact Hll regions in the Milky Way 


(Walsh et al. 1999). As an independent check, the H 2 density calculated by Turner et al. 
(2015) corresponds to p gas ~1.56xl0 -19 g cm -3 , or a dust/gas ratio~0.0045 for the cloud of 


cluster 11, which is only slightly lower than the values measured for nearby galaxies (Draine 


et al. 2007). This ratio could be lower still, since the H 2 density is derived from the CO(3-2) 
transition, which yields a lower limit to the actual molecular gas density. High-resolution 
HI maps of the region only show HI absorption, which sets a lower limit of 5 x 10 2 ° cm 


20 2 tQ 


the HI column density (Kobulnicky & Skillman 1995) 


The virial mass contained in a region about 43 pcx23 pc, roughly the size of the region 
displayed in either panel of Figure [ 2 J is about 1.2-1.3xl0 6 M 0 (Turner et al. 2015, scaled 


to our distance). If clusters 5 and 11 provide most of the stellar mass in the region, the star 


formation efficiency is SFE~0.25-0.30, about a factor 2 lower than the estimate 13 of Turner 


et al. (2015). The SFE we derive is comparable to the SFE of clusters in our own Milky 
Way, but larger than the SFE~5%-10% measured when including the entire molecular cloud 


(Lada & Lada 2003). Whether or not clusters 5 and 11 will remain bound depends on both 


the value of the local (cluster-size) SFE and the ratio of the gas removal to the crossing 


timescales (Parmentier et al. 2009) 


If they remain bound and do not lose significant mass, clusters 5 and 11 are massive 
enough that they could be progenitors of globular clusters. Globular clusters, like those in the 
halo of our own Milky Way and other galaxies, have masses in the range ?al0 4 -10 6 M 0 (Fall, 
Chandar & Whitmore 2009]). The mass loss of an evolving cluster depends on a number of 
factors, including the nature/origin of the second population (Schaerer & Charbonnel 2011), 
but under most scenarios, both clusters 5 and 11 would retain sufficient mass to remain 
within the range of globular clusters. Clusters 5 and 11 are two of the youngest among 
the very massive clusters (super-star-clusters) detected in nearby galaxies, such as those in 


He2- 

10 ( 

Kobulnicky & Johnson 1999), NGC1569 ( 

Hunter et al. 

2000 

), M82 ( 

Smith et al. 

2006 

), SBS0335-052 ( 

Johnson, Hunt & Reines 

2009 

), the Antennae ( 

Whitmore et al. 

2010 

), 


and NGC1705 (Martins et al. 2012), to name a few. As such, they provide important case 


studies to test theories of multiple populations in globular clusters (e.g. de Mink et al. 2009 


Bastian, Cabrera-Ziri, & Salaris 2015). 


13 The star formation efficiency decreases to SFE~0.15, if we adopt a top-heavy stellar IMF with a lower 
cut-off of 1 M 0 ; this cut-off is lower than the 3 M 0 proposed by 


Turner et al. 


(2015). 
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6.6. The Starburst in NGC 5253 


The youth of the two clusters in the radio nebula has an interesting implication for the 
estimate of the SFR in this galaxy from its TIR luminosity, L T j f{ =?).7x 10 42 erg s _1 . The 
majority of the TIR emission in NGC 5253 originates from the same region, as suggested by 
the point-like emission in the Spitzer images (Dale et al. 2009) and by mid-IR imaging (Go: 


jian, Turner & Beck 2001). If we use a scaling factor between L(TIR) and SFR appropriate 
for a population as young as 2 Myr (Calzetti 2013), we obtain SFR(TIR)=0.31 M 0 yr -1 , 
a factor almost 3 larger than what would be inferred using the standard calibration. The 
higher SFR(TIR) is consistent with the range of SFRs derived from the free-free emission, 
SFR~0.3-0.36 M 0 yr^ 1 , suggesting that, within the uncertainties, we see no obvious indica¬ 
tion for significant direct absorption of ionizing photons by dust. Adding the UV, which orig¬ 
inates in separate regions from the IR, we get a total SFR for NGC 5253 SFR=0.37 M 0 yr -1 , 
when using a timescale of 10 Myr for the SFR(UV), which is appropriate for the youth of 
the UV-emitting area. 

Moving away from the central radio nebula, the clusters increase in age for increasing 
distance from clusters 5 and 11 (Figures [TT7| and [IT]) . Clusters 1-4 have best-fit ages around 
5-6 Myr, while clusters 7-10 are consistent with ages ~9-15 Myr. The ‘association’ of 
clusters 1-4, all with comparable ages and spread over a small area of ~170 pc 2 , agrees with 


the micro-level hierarchy picture discussed by Efremov & Elmcgreen (1998); these authors 


find that in the Magellanic Clouds small regions form stars over a short period of time. UV- 


bright star clusters to the North of the radio nebula have ages in the range 4-8 Myr (Tremonti 


et al. 2001). de Grijs et al. (2013) also recovers typically young ages for the star clusters 
in the central region of NGC 5253, between a few 10 6 yr and a few 10' yr. This suggests a 
progressive concentration of the most recent event of star formation from a larger, ~300 pc, 
region to the smaller area, ~20 pc, of the radio nebula. Covering 300 pc in ~15 Myr requires 
a propagation velocity of the perturbation of about 20 km s -1 . This is larger than the typical 
sound speed of a few krn/s of the ISM, but is consistent with the velocities of a few tens 
of krn/s expected for the propagation of turbulence in a multi-phase ISM (e.g. 


Bournaud 


et al. 2011). External triggers, such as the recent-past interaction with M83 and/or gas 


infall, are required for such large-scale propagating disturbances, although internal triggers 


(e.g., shocks from earlier events that produce sequential star formation, see, Whitmore et al. 


2010) may also play a role. The region within the central ^300 pc of NGC 5253 is likely to 


have been the stage of one continuous and connected event of star formation over at least the 
past ~15 Myr. This is in agreement with other estimates that indicate relatively constant 


levels of SFR over the past few hundred Myr (McQuinn et al. 2010 Harbeck, Gallagher & 


Crnojevic|[2012), although these studies cannot discriminate between continuous or sporadic 


star formation over this timescale. 
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The energetics that drive the ISM kinematics in NGC5253 are also consistent with a 
continuous star formation event over at least this timescale. The measured Ha luminosity of 
the shocked gas, L(H a) s h oc k ~3xl0 39 erg s -1 (Hong et al. 2013, , rescaled to our distance), 


corresponds to a mechanical luminosity of L mec h ~7xl0 4 ° erg s 1 ( 

dich et al. 2010), for 

radiative shocks and the shock velocity measured by 

Marlowe et al. 

(1995 

). For constant 


star formation and our adopted Kroupa IMF, the mechanical energy requirements can be 
satisfied with SFR~O.2M 0 yr -1 over the past 10-15 Myr. The five star clusters in this age 
range, 6-to-10, correspond to a cumulative mass of ~1.6xl0 5 M 0 (about 1/2 of the sum of 
the masses of 5 and 11), and may have contributed, in the past, about 1/3 of the mechanical 
energy needed to support the shocks. The rest of the energy requirement has come from 
the remaining young star clusters and diffuse stars in the region. When they will have aged 
enough to start producing supernova explosions, clusters 5 and 11 will be able to supply 
about 2/3 of the mechanical energy requirement of the present-day shocks. 


7. Summary and Conclusions 


The combination of new and archival HST imaging data of the central 300 pc in the 
starburst galaxy NGC5253 has enabled us to derive, with unprecedented accuracy, ages, 
masses, and extinctions for the brightest among the star clusters in the starburst. The HST 
data cover the full wavelength range from ~1500 A to 1.9 /jrn in 10 continuum bands and 
3 narrow bands, which include the emission lines of Ha, P/3, and Pa. The 11 star clusters 
analyzed here include the two young clusters located within the central dusty radio nebula, 
in addition to nine clusters distributed across the UV-bright starburst. 


The multi-wavelength SED of each cluster has been fit with models that include stellar 
populations, gas emission, and dust attenuation (and, for cluster 11, dust emission). Models 
that include only foreground dust are sufficient to produce excellent fits for the 10 bright 
clusters, while the heavily dust extincted cluster 11 requires a geometry that includes a 
combination of foreground and mixed dust. 


Because of its location within a region of high optical depth, this is the first time the 
SED of cluster 11 has been fit shortward of the J-band. This is also the first time extensive 
stellar continuum photometry, from the UV to the nearIR, is available for a natal cluster. 
Albeit counter-intuitive, cluster 11 is detectable in all optical and UV bands, mainly because 
dust mixed with stars dims the light more readily than reddening it (Calzetti 2001). Like 
for all other clusters analyzed in this paper, the SED of cluster 11 can be fit with a standard 


Kroupa (2001) IMF in the mass range 0.1-120 M 0 , and does not require truncations or other 


modifications. The resulting best-fit age, mass, and total dust attenuation of cluster 11 yield 
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a self-consistent picture of this source, also in agreement with independent constraints. 


All clusters studied here are younger than ~15 Myr, in agreement with previous results, 
with a systematic trend for the age to become younger when moving from the outskirts of 
the starburst towards its center, where the radio nebula is located. The picture that emerges 
is for the current burst of star formation to have been a continuous event over the past 
~15 Myr, and to have concentrated inward, forming a hierarchical ‘age’ structure similar 


to those discussed by Efremov & Elmcgreen (1998) for the LMC. Dwarf starburst galaxies, 


including NGC 5253, do tend, in fact, to display a single dominant region of connected star 
formation hierarchical in structure (Elmegreen et al. 2014). SFR indicators calibrated for 


the youth of the starburst yield SFR~O.4 0 yr 1 for the galaxy, or about three times higher 
than using standard SFR calibrations. 


The nine star clusters located outside the radio nebula have masses in the range 0.5- 
5xl0 4 M 0 , and dust attenuations Ay <1 mag. The relatively low Ay value is consistent 
with the location of the clusters in the UV-bright region of the starburst. Other indicators 
show dust attenuation to be low in this region in general. The nine clusters provide today a 
minimal contribution, less than 2%, to the total ionizing photons in the galaxy. 

The two most massive clusters in our sample, clusters 5 and 11, are both located inside 
the dusty radio nebula. Their large masses, 7.5xlO 4 M 0 and 2.5xlO 5 M 0 , respectively, are 
consistent with having been born in a high-density environment. Both clusters are extremely 
young, with best fit ages 1±1 Myr. Neither has reached a stage when supernovae are starting 
to affect the local environment, including blowing out the natal dust cloud, which accounts 
for their significant dust attenuation. Energy balance between the ionizing photon rate 
recovered from the free-free emission and the total far-infrared luminosity suggests that, 
albeit buried in a significant dust cloud, the clusters in the radio nebula do not suffer from 
significant direct absorption of ionizing photons by dust. 

Both clusters are behind a dust layer of Ay <2 mag, which is likely the outer dust 
shell of the radio nebula. Cluster 5 is behind only this dust layer, which is still sufficiently 
transparent to enable this cluster to be the Ha peak in the galaxy. Cluster 11 is attenuated 
by an additional dust cloud with front-to-back optical depth Ay ~49 mag, mixed with 
its stars. This dust cocoon has many characteristics typical of those observed in Milky 
Way Ultracompact HII regions, including a maximum temperature of ~ 1,100 K for the dust 
closest to the massive stars. It is also responsible for about 60%, 65%, and 90% of the 
emission observed in the J, H, and K bands, respectively. 

The mass of cluster 11 places it a factor ^2-4 below the ~0.6-1.2xl0 6 M 0 value inferred 
in previous studies. Except for its ‘dusty’ condition, owed presumably to its youth, this 
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cluster’s other characteristics are relatively normal. For instance, it is within 2-3a of the 
relation between the most massive clusters and SFR in galaxies, and is similar to massive 
clusters recovered in other dwarf galaxies. Its mass is in-between those of the scaled OB 


association NGC604 in M33 ( 

Maiz-Apcllaniz 

2001 

) and the massive cluster NGC346 in the 

SMC ( 

Portegies Zwart et al. 

2010 

). For the level of dust attenuation and for its youth, 


cluster 11 is reminiscent of the super-star-cluster WS 80 in the Antennae galaxy (Whitmore 


& Zhang 2002), although it is also about 16 times less massive. Based on the self-consistency 


of several indicators, we do not believe to be missing a large fraction of the mass of cluster 11 . 


Clusters 5 and 11 provide about half of the ionizing photons in the radio nebula, and 
about one third of the ionizing flux in the entire galaxy. The remaining 50% of the nebula’s 
ionizing photons are diffuse, suggesting leakage of ionizing photons from the immediate sur¬ 
roundings of the two clusters into the nebula’s region. Once they start producing supernovae, 
clusters 5 and 11 will supply about 2/3 of its mechanical luminosity requirements to support 
the current level of shocks. Thus, by themselves, the two clusters can sustain a significant 
fraction of the energy requirements of NGC 5253. The remaining energy needs to be supplied 
by other sources, which may include: other stars and star clusters within the area of the 


radio nebula, Very Massive Stars (< 

Lrowther et al. 2010 

Doran et al. 2013 

), and/or rotating 

massive stars ( 

Leitherer et al. 

201 4 

). The potential presence of VMSs is consistent with the 


detection of Wolf-Rayet features of WN type in the region of the radio nebula. Follow-up 
high-spatial-resolution spectroscopy of these clusters has the potential to reveal signatures 
of rare star formation products and address this open question. 


Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at 
the Space Telescope Science Institute, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under NASA contract NAS 5-26555. These observations 
are associated with program 7 ^ 13364. Support for program 7 ^ 13364 was provided by NASA 
through a grant from the Space Telescope Science Institute. 

Based also on observations made with the NASA/ESA Hubble Space Telescope, and 
obtained from the Hubble Legacy Archive, which is a collaboration between the Space Tele¬ 
scope Science Institute (STScI/NASA), the Space Telescope European Coordinating Facility 
(ST-ECF/ESA) and the Canadian Astronomy Data Centre (CADC/NRC/CSA). 

This research has made use of the NASA/IPAC Extragalactic Database (NED) which 
is operated by the Jet Propulsion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Administration. 
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tinguished Visitor at the Institute of Astronomy, University of Cambridge (UK), and an 
















39 


Overseas Fellow at the Churchill College (Cambridge, UK). D.C. acknowledges the kind hos¬ 
pitality of both the Institute and the College. A.S.E. was supported by the Taiwan, R.O.C. 
Ministry of Science and Technology grant MoST 102-2119-M-001-MY3. M.F. acknowledges 
support by the Science and Technology Facilities Council [grant number ST/L00075X/1], 
D.A.G. kindly acknowledges financial support by the German Research Foundation (DFG) 
through grant GO 1659/3-2. E.Z. acknowledges research funding from the Swedish Research 
Council (project 2011-5349). 


REFERENCES 

Alonso-Herrero, A., Takagi, T., Baker, A.J., Rieke, G.H., Ricke, M.J., Imanishi, M., & 
Scoville, N.Z. 2004, ApJ, 612, 222 

Anders, P., Kotulla, R., de Grijs, R., & Wicker, J. 2013, ApJ, 778, 138 

Asplund, M., Grevesse, N., Sauval, A.J., & Scott, P. 2009, ARA&A, 47, 481 

Bastian, N., Cabrera-Ziri, I., & Safaris, M. 2015, MNRAS, 449, 3333 

Beck, S.C., Turner, J.L., Ho, P.T.P., Lacy, J.H., & Kelly, D.M. 1996, ApJ, 457, 610 

Billett, O.H., Hunter, D.A., & Elmegreen, B.G. 2002, AJ, 123, 1454 

Bohlin R.C., Savage B.D., Drake J.F. 1978, ApJ, 224, 132 

Bournaud, F., Chapon, D., Teyssier, R., Powell, L.C., Elmegreen, B.G., Elmegreen, D.M.., 
Due, P.A., Contini, T., Epinat, B., & Shapiro, K.L. 2011, ApJ, 730, 4 

Bresolin, F. 2011, ApJ, 729, 56 

Caldwell, N., & Phillips, M.M. 1989, ApJ, 338, 789 

Calzetti, D. 2001, PASP, 113, 1449 

Calzetti, D. 2013, in ‘Secular Evolution of Galaxies’, J.Falcon-Barroso & J.H. Knapen eds. 
(Cambridge, UK: Cambridge LIniversity Press), 419 

Calzetti, D., Kinney, A.L., & Storchi-Bergmann, T. 1994, ApJ, 429, 582 

Calzetti, D., Meurer, G.R., Bohlin, R.C., Garnett, D.R., Kinney, A.L., Leitherer, C., & 
Storchi-Bergmann, T. 1997, AJ, 114, 1834 

Calzetti, D., Conselice, C.J., Gallagher, J.S. III, & Kinney, A.L. 1999, AJ, 118, 797 



40 


Calzetti, D., Armus, L., Bohlin, R.C., Kinney, A.L., Koornneef, J., & Storchi-Bergmann, T. 
2000, ApJ, 533, 682 

Calzetti, D., Harris, J., Gallagher, J.S., Smith, D., Conselice, C.J., Homeier, N., & Kewley, 
L. 2004, AJ, 127, 1405 

Calzetti, D., Lee, J.C., Sabbi, E., Adamo, A., Smith, L.J., Andrews, J.E., Ubeda, L., Bright, 
S.N., Thilker, D., Aloisi, A., et al. 2015, AJ, 149, 51 

Campbell A., Terlevich R., & Melnick J., 1986, MNRAS, 223, 811 

Chandar, R., Leitherer, C., Tremonti, C.A., Calzetti, D., Aloisi, A., Meurer, G.R., & de 
Mello, D. 2005, ApJ, 628, 210 

Cresci G., Vanzi L., & Sauvage M. 2005, A&A, 433, 447 

Cervino, M., & Luridiana, V. 2004, A&A, 413, 145 

Cook, D.O., Dale, D.A., Johnson, B.D., van Zee, L., Lee, J.C., Kennicutt, R.C., Calzetti, 
D., Staudaher, S.M., & Engelbracht, C.W. 2014, MNRAS, 445, 899 

Crowther, P.A., Schnurr, O., Hirschi, R., Yusof, N., Parker, R.J., Goodwin, S.P., & Kassim, 
H.A. 2010, MNRAS, 408, 731 

Dale, D.A., Cohen, S.A., Johnson, L.C., Schuster, M.D., Calzetti, D., Engelbracht, C.W., 
Gil de Paz, A., Kennicutt, R.C., Lee, J.C., Begum, A., et al. 2009, ApJ, 703, 517 

Davidge, T.J. 2007, AJ, 134, 1799 

de Grijs, R., Anders, P., Zackrisson, E., & Ostlin, G. 2013, MNRAS, 431, 2917 

de Mink, S.E., Pols, O.R., Langer, N., & Izzard, R.G. 2009, A&A, 507, LI 

Doran, E.I., Crowther, P.A., de Koter, A., Evans, C.J., McEvoy, C., Walborn, N.R., Bastian, 
N., Bestenlehner, J.M., Graefener, G., Herrero, A., et al. 2013, A&A, 558, A134 

Draine, B.T. 2003, ARA&A, 41, 241 

Draine, B.T. 2004, in ‘The Cold LIniverse’, eds. A.W. Blain, F. Combes, B.T. Draine, D. 
Pfenniger & Y. Revaz, Saas-Fee Advanced Course 32 (Springer-Verlag), 213 

Draine, B.T., Dale, D.A., Bendo, G., Gordon, K.D., Smith, J.D.T., Armus, L., Engelbracht, 
C.W., Hclou, G., Kennicutt, R.C., Li, A., et al. 2007, ApJ, 663, 866 

de Mink, S.E., Sana, H., Lange, N., Izzard, R.G., & Schneider, F.R.N. 2014, ApJ, 782, 7 



-41 - 

Eldridge, J.J., & Stanway, E.R. 2009, MNRAS, 400, 1019 
Efremov, Y.N., & Elmegreen, B.G. 1998, MNRAS, 299, 588 

Elmegreen, D.M., Elmegreen, B.G., Adamo, A., Aloisi. A. Andrews, J., Annibali, F., Bright, 
S.N.,, Calzetti, D., Cignoni, M., Evans, A.S., et al. 2014, ApJ, 787, L15 

Fall, S.M., Chandar, R., & Whitmore, B.C. 2009, ApJ, 704, 453 

Ferland, G.J., Porter, R.L., van Hoof, P.A.M. , Williams, R.J.R. , Abel, N.P. , Lykins, M.L. 
, Shaw, G., Henney, W.J., & Stancil, P.C. 2013, RMexAA, 49, 1 (The 2013 Release 
of Cloudy) 

Fitzpatrick, E.L. 1999, PASP, 111, 63 

Freedman, W.L., Madore, B.F., Gibson, B.K., Ferrarese, L., Kelson, D.D., et al. 2001, ApJ, 
553, 47 

Gazak, J.Z., Bastian, N., Kndritzki, R.-P., Adamo, A., Davies, B., Plez, B., & Urbaneja, 
M.A., 2013, MNRAS, 430, L35 

Girardi, L., Bressan, A., Bertelli, G., & Chiosi, C. 2000, A&ASnpp., 141, 371 
Gorjian, V., Turner, J.L., & Beck, S.C. 2001, ApJ, 554, L29 

Grimes, J.P., Heckman, T., Aloisi, A., Calzetti, D., Leitherer, C., Martin, C.L., Meurer, G., 
Sembach, K., & Strickland, D. 2009, ApJS, 181, 272 

Harbeck, D., Gallagher, J.S. Ill, & Crnojevic, D. 2012, MNRAS, 422, 629 

Harris, J., Calzetti, D., Gallagher, J.S., Smith, D.A., & Conselice, C.J. 2004, ApJ, 603, 503 

Hong, S., Calzetti, D., Gallagher, J.S., III; Martin, C.L., Conselice, C.J., & Pellerin, A. 2013, 
ApJ, 777, 63 

Hunter, D.A., O’Connell, R.W., Gallagher, J.S., & Smecker-Hane, T.A. 2000, AJ, 120, 2383 
Johnson, K.E., Vacca, W.D., Leitherer, C., Conti, P.S., & Lipscy, S.J. 1999, AJ, 117, 1708 
Johnson, K.E., Indebetouw, R., Watson, C., & Kobulnicky, H.A. 2004, AJ, 128, 610 
Johson, K.E., Hunt, L.K., & Reines, A.E. 2009, AJ, 137, 3788 
Karachentsev I. D. et ah, 2007, AJ, 133, 504 



-42- 

Kauffmann, J., Filial, T., & Zhang, Q. 2013, A&A, 765, L35 

Kennicutt, R.C., & Evans, N.J. 2012, ARA&A, 50, 531 

Kennicutt, R.C., Lee, J.C., Funes, S.J. , Sakai, S., & Akiyama, S. 2008, ApJS, 178, 247 

Kobayashi, H., Kimura, H., Watanabe, S., Yamamoto, T., & Muller, S. 2011, Earth, Planets, 
and Space, 63, 1067 

Kobulnicky H.A., & Skillman, E.D., 1995, ApJ, 454, L121 

Kobulnicky H.A., Skillman E.D., Roy J.-R.,Walsh J.R., & Rosa M.R., 1997, ApJ, 277, 679 

Kobulnicky H.A., & Johnson, K.E., 1999, ApJ, 527, 154 

Koehler, K., Langer, N., de Koter, A,; de Mink, S.E., Crowther, P.A., Evans, C.J., Grafener, 
G., Sana, H., Sanyal, D., Schneider, F.R.N., & Vink, J.S. 2015, A&A, 573, A71 

Koekemoer, A.M., McLean, B., McMaster, M., & Jenkner, H. 2006, in ‘The 2005 HST 
calibration workshop: Hubble after the transition to two-gyro mode’, Proceedings of 
a workshop held at the Space Telescope Science Institute, Baltimore, Maryland; A. 
M. Koekemoer, P. Goudfrooij, and L. L. Dressel eds., p. 384 

Kreckel, K.,Groves, B., Schinnerer, E., Johnson, B.D., Aniano, G., Calzetti, D., Croxall, 
K.V., Draine, B.T., Gordon, K.D., Crocker, A.F., et al. 2013, ApJ, 771, 62 

Kroupa, P. 2001, MNRAS, 322, 231 

Lada, C.J., & Lada, E.A. 2003, ARA&A, 41, 67 

Lee, M.G., Chandar, R., & Whitmore, B.C. 2005, AJ, 130, 2128 

Leitet, E., Bergvall, N., Hayes, M., Linne, S., & Zackrisson, E. 2013, A&A, 553, A106 

Leitherer, C., Schaerer, D., Goldader, J.D., Gonzalez Delgado, R.M, Robert, C., Kune, D.F., 
de Mello, D.F., Devost, D., & Heckman, T.M. 1999, ApJS, 123, 3 

Leitherer, C., Eckstrom, S., Meynet, G., Schaerer, D., Agienko, K.B., & Levesque, E.M., 
2014, ApJS, 212, 14 

Lopez-Sanchez, A.R., Esteban C., Garcia-Rojas, J., Peimbert, M., & Rodriguez, M., 2007, 
ApJ656, 168 

Lopez-Sanchez, A.R., Koribalski, B.S., van Eymeren J., Esteban C., Kirby E., Jerjen H., & 
Lonsdale N., 2012, MNRAS, 419, 1051 



Maiz-Apellaniz, J. 2001, ApJ, 563, 151 

Maoz, D., Barth, A.J., Sternberg, A., Filippenko, A.V., Ho, L.C., Macchetto, F.D,; Rix, 
H.-W., & Schneider, D.P. 1996, AJ, 111, 2248 

Marlowe A.T., Heckman T.M., Wyse R.F.G., & Schommcr R., 1995, ApJ, 438, 563 

Martin, C.L., 1998, ApJ, 506, 222 

Martins, F., Foerster Schreiber, N.M., Eisenhauer, F., & Lutz, D. 2012, A&A, 547, A17 
Martin-Hernandez, N.L., Schaerer, D., & Sauvage, M. 2005, A&A, 449, 2005 
McQuinn, K.B.W. et ah, 2010, ApJ, 724, 49 
Meier, D.S., Turner, J.L., & Beck, S.C. 2002, AJ, 124, 877 

Meurer, G.R., Heckman, T.M., Leitherer, C., Kinney, A., Robert, C., & Garnett, D.R., 1995, 
AJ, 110, 2665 

Meynet, G., Schaller, G., Schaerer, D., & Charbonnel, C., 1994, A&ASupp., 103, 97 
Monreal-Ibero A., Vlchez J.M., Walsh J.R., & Munoz-Tunon C., 2010, A&A, 517, 27 
Monreal-Ibero, A., Walsh, J. R., & Vlchez, J. M. 2012, A&A, 544, A60 
Monreal-Ibero, A., Walsh, J.R., Westmoquette, M.S., & Vlchez, J.M. 2013, A&A, 553, A57 
Moorwood, A.F.M., & Glass, I.M. 1982, A&A, 115, 84 
Moustakas, J.,& Kennicutt, R.C., 2006, ApJS, 164, 81 

Parmentier, G., Goodwin, S.P., Kroupa, P., & Baumgardt, H. 2009, ApJ, 678, 347 
Pcllcrin, A., & Robert, C. 2007, MNRAS, 381, 288 
Popescu, B., & Hanson, M.M. 2010, ApJ, 713, L21 

Portegies Zwart, S.F., McMillan, S.L.W., & Gieles, M. 2010, ARA&A, 48, 431 
Reines, A.E., Nidever, D.L., Whelan, D.G., & Johnson, K.E. 2010, ApJ, 708, 26 
Rich, J.A., Dopita, M.A., Kewley, L.J., & Rupke, D.S.N. 2010, ApJ, 721, 505 
Rieke, G.H., Lebofsky, M.J., & Walker, C.E. 1988, ApJ, 325, 679 
Sakai, S., Ferrarese, L., Kennicutt, R.C., & Saha, A. 2004, ApJ, 608, 42 



- 44 - 

Sana, H., de Mink, S.E., de Koter, A., Langer, N., Evans, C.J., Gicles, M., Gosset, E., Izzard, 
R.G., Le Bouquin, J.-B., & Schneider, F.R.N., 2012, Science, 337, 444 

Schaerer ,D., Contini, T., Kunth ,D., & Meynet, G., 1997, ApJ, 481, L75 

Schaerer, D., & Charbonnel, C. 2011, MNRAS, 413, 2297 

Schlafly, E.F., & Finkbeiner, D.P., 2011, ApJ, 737, 103 

Schneider, F.R.N., Izzard, R.G., de Mink, S.E., Langer, N., Stolte, A., de Koter, A., Gvara- 
madze, V.V., Hussmann, B., Liermann, A., & Sana, H. 2014, ApJ, 780 117 

Schneider, F.R.N., Izzard, R.G., Langer, N., & de Mink, S.E. 2015, ApJ, 805, 20 

Smith, L.J., Westmoqnette, M.S., Gallagher, J.S., O’Connell, R.W., Rosario, D.J., & de 
Grijs, R, 2006, MNRAS, 370, 513 

Storchi-Bergmann, T., Kinney, A. L., & Challis, P. 1995, ApJS, 98, 103 
Strickland, D.K., & Stevens, I.R. 1999, MNRAS, 306, 43 

Summers, L.K., Stevens, I.R., Strickland, D.K., & Heckman, T.M. 2004, MNRAS, 351, 1 

Thim F., Tammann G.A., Saha A., Dolphin A., Sandage A., Tolstoy E., & Labhardt L., 
2003, ApJ, 590, 256 

Turner, J.L., & Welch, W.J. 1984, ApJ, 287, L81 
Turner, J.L., Beck, S.C., & Hurt, L.R. 1997, ApJ, 474, Lll 
Turner, J.L., Beck, S.C., & Ho, P.T.P. 2000, ApJ, 532, L109 
Turner, J.L., & Beck, S.C. 2004, ApJ, 602, L85 

Turner, J.L., Beck, S.C., Bendford, D.J., Consiglio, S.M., Ho, P.T.P., Kovacs, A., Meier, 
D.S., & Zhao, J.-H. 2015, Nature, 519, 331 

Tremonti, C.A., Calzetti, D., Leitherer, C., & Heckman, T.M. 2001, ApJ, 555, 322 

van den Bergh, S. 1980, PASP, 92, 122 

Vanzi, L., & Sauvage, M. 2004, A&A, 415, 509 

Vazquez, G.A., & Leitherer, C. 2005, ApJ, 621, 695 

Walsh, J. R., & Roy, J.-R. 1989, MNRAS, 239, 297 



Walsh, A.J., Burton, M.G., Hyland, A.R., & Robinson, G. 1999, MNRAS, 309, 905 
Westmoquette, M.S., James, B., Monreal-Ibero, A., & Walsh, J.R. 2013, A&A, 550, A88 
Whitmore, B.C., & Zhang, Q. 2002, AJ, 124, 1418 

Whitmore, B.C., Chandar, R., Schweizer, F., Rothberg, B., Leitherer, C.,Rieke, M., Rieke, 
G., Blair, W.P., Mengel, S., & Alonso-Herrero, A., 2010, AJ, 140, 75 

Whitmore, B.C., Chandar, R., Kim, H., Kaleida, C., Mutchler, M., Stankiewicz, M., Calzetti, 
D., Saha, A., O’Connell, R., Balick, B., et ah 2011, ApJ, 729, 78 

Whitmore, B.C., Brogan, C., Chandar, R., Evans, A., Hibbard, J., Johnson, K., Leroy, A., 
Privon, G., Remijan, A., & Sheth, K., 2014,a ApJ, 795, 156 

Whitmore, B.C., Chandar, R., Bowers, A.S., Larsen, S., Lindsay, K., Ansari, A., & Evans, 
J. 2014b, AJ, 147, 78 

Yusof, N., Hirschi, R., Meynet, G., Crowther, P.A., Ekstrom, S., Frischknecht, U., Georgy, 
C., Abu Kassim, H., & Schnurr, O. 2013, MNRAS, 433, 1114 

Zackrisson, E., Bergvall, N., Oloffson, K., & Siebert A. 2001, A&A, 375, 814 

Zackrisson, E., Rydberg, C.-E., Schaerer, D., Ostlin, G., & Tuli, M. 2011, ApJ, 740, 13 

Zastrow, J., Oey, M.S., Veilleux, S., McDonald, M., & Martin, C.L. 2011, ApJ, 741, L17 

Zastrow, J., Oey, M.S., Veilleux, S., & McDonald, M. 2013, ApJ, 779, 76 


This preprint was prepared with the AAS IATgX macros v5.2. 



CO 

10 

CN 

iO 

O 

O 

Z 


C/3 

CD 


b€ 

cb 


a 

HH 


EH 

co 


K 


CD 

4b 

h-h 


03 

O 


S-i 

CD 

-4> 

o 

Ctf 

S-h 

<3 

O 


43 

Eh 


bjO 

O 

f—i 

Ph 


'O # 
03 
X? 

O 

CD 

"§ 

P 



£ 

CD 

> 

<-M 

o 

2 

2 

Ph 


X ^ 


0 ) 

Ph . 

bb ' 
d 


c6 

O 

N 


CD 

X 

P 


o 


d 

CD 


05 



05 

05 

05 

05 

05 

CO 

CO 





b- 

CO 

CO 

o 

o 

O 

O 

O 

o 

O 

05 

05 

05 

05 

kO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CM 

CM 

t-H 

rH 

rH 

t-H 

i—H 

CO 

CO 

o 

o 

o 

o 

o 

CM 

CM 

CM 

CM 

CM 

CM 

t—1 

1—1 

rH 

rH 

1 

rH 

T—1 

rH 

I - 

^6 

1 t-H 

i r 

b- 

1 

b- 

b- 

b- 

1 

O 

6 

6 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

c 

a 

o 

o 

o 

a 

a 

o 

a 

a 

a 

0 

a 

o 

o 

!- 

00 

00 

o 

o 

o 

o 

o 

co 

CO 





O 

CM 

1 

CM 

1 

CM 

CH 

1 

1 

CN 

1 

CM 

l 

CM 

CN 

1 

o 

o 

Q 

o 

CO 

00 

00 

CM 

CM 

CM 

CM 

CM 

Ib- 

b!- 

T—1 

t—H 

t-H 

rH 

o 

o 

o 

o 

O 

O 

o 

O 

CD 

o 

o 

O 

o 

o 

05 

CO 

CO 

CO 

co 

CO 

CO 

CO 

1 

t-H 

i 

rH 

00 

00 

00 

00 

o 

t-H 

t—H 

o 

o 

O 

o 

o 

t-H 

t-H 

05 

05 

05 

05 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

05 

05 

05 

05 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

t—H 

t-H 

t-H 

t—H 

^__ 

_ 

_ 

_ 

^^ 

__ N 

_ 

_ 

_ 

^^ 

_ 

_ 

_ 

^^ 


CO 

CO 






CO 

CO 





X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

kO 

CO 

CM 

05 

C5 

o 

CD 

CM 

CM 

CM 





CO 

t-H 

00 


kO 

o 

CO 

05 

05 

05 

CM 

CM 

co 

CO 

CO 

00 

b- 



CM 

x — 


05 

05 

^ 

x —' 

"■— 

s — 


s —■" 



S ^' 

s — 

o 

CD 

T—H 


CD 

C5 

CD 

CD 

o 

o 

C5 

CD 

CD 

o 

o 

00 

CO 

CO 

CO 

CO 

cd 

cd 

o 

00 

CO 

CO 

O 

CD 


CO 

05 

05 

LO 

kO 

CO 


Tt 

05 

o 

O 

CM 

CO 

t— 

fP 



CM 

CM 

CO 


CO 

b- 

CO 

00 



05 

05 





CM 

CM 

CM 

T—1 





kO 















t-H 





CM 

r —H 

00 

b- 

CD 

00 

b— 

2 

o 

CD 

O 

C5 

00 

CD 

00 

C5 


CM 

t-H 

05 

CO 

kO 


CM 

CM 

t-H 

b4 

cd 

CO 

rH 

kO 

CO 

t-H 

b- 

00 

T—i 

CO 

CO 

05 

b- 


00 


b- 

CO 

CO 

CO 

kO 

kO 

t-H 

LO 

00 

CM 

o 

b- 

05 

l—H 

CM 

CO 

CO 

Hf 

kO 

CO 

00 

t-H 

CM 

t—H 

CO 

00 

00 









t-H 

t-H 

t—H 

t-H 

t-H 

t-H 





















HH ■ 














HH 




























+ 



^- s 



^— s 











'Q. 



a 


Oh 

DP 









d, 



Ph, 


^ ^ ^ ^ S 

& 

^ ^ ^ ^ ^ 

2 ; 

£ 

kO 

kO 

CO 

o 

kO 

o 

00 

TJH 

o 

00 

o 

o 


o 

CM 

15- 

CO 

CO 

CO 

kO 

kO 

t-H 

t-H 

CM 

t—H 

CO 

oo 

05 

rH 

CM 

CO 

CO 


kO 

CO 

00 

t-H 

t-H 

t-H 

rH 

t—H 

t-H 

Ph 

Ph 

Ph 

Ph 

Ph 

Ph 

Ph 

Ph 

Ph 


Ph 

Ph 

Ph 



CM 







co 







CO 







CO 






t—H 

rH 


co 





t-H 


05 




CO 

X 


CM 





X 


t—H 




X 

CM 


X 





CO 


X 




kO 

CO 


05 





CM 


05 




CO 

i-H 


CM 





t-H 


t—H 




b- 

o 


o 





CM 


CO 




05 

b- 


CO 





CM 


t—H 




o 

o 


CD 





CD 


o 




o 

o 


CD 











o 

o 


LO 





O 


kO 




CO 



CM 





CO 


b- 




o 

CD 


CD 





t-H 


CD 




o 

o 


CD 





CD 


CD 














CM 














0 





m 









HH 




o 

Pb 

m 

HH 

> 

h3 

CO 


o 

pp 





PP 

HH 

CO 


m 

O 

2 

o 

HH 




m P 
O 4 

<c 4 


m 

O 

c 





o 

Ph 

4 






a m 

> ^ 

s 

cp 
CO Ph 
Ph <u 
^ 43 

g h 


05 

4 

o 

05 

> 


CD 

d 

s5 

^ _ 

c3 CG 
05 rt 
4 m 
Oh 05 
CO 4b 


O 

Oh 

05 

43 


r Tb 

05 


d 

O 

43 

-4 

05 

2 

CD 


05 

X 


w Cl 

C/3 ^ 

05 

n 5b 

Oh 05 

X * 
0 > 2 
s s 

O XJ 




4b 

bjO 


d ^ 
bJO -4 


d 

d 

05 

4b 


CD 

4 

r rD 

O 

4 

Oh 


05 05 

rb bjO 
o 2 


05 

X 

‘Bh 

05 

> 


d bb 

C/3 

+3 05 

4b 
4b 4= 

~4 „ 

•d C/3 
S 05 
C/3 

4 o3 
05 CD 
4b 
-4 
05 
bjO 
O 

-4 

d~ 

.2 

’-4 

o3 
b 


a 

05 

05 

4b 


d 

o 


CD 

c3 


4 

o 


d 

d 

.3 

*-4 

d 

o 

CD 


-4 <3 

05 r-H 

bo 05 
H -4 
CO C/3 

-* <u 

£ 3 

b~ 4 

00 o 

t-H 4h 

Ph 

05 

C/3 

c 3 

as .a 


05 

bO 

cO 


00 
CM 

t-H 

fc T3 

^ s 

„ cS 

£ 4 
00 1 
LO t> 

CO o 

s 

a s 

<u g 

•3 ° 

c£ Oi 

T—I 

"O Ph 

a <u 

H2 H 

o 4 

4 <p 

> 


cO 

d 

05 

4b 

Eh 


CD 

05 

Oh 


°<1 
I O 
, LO 
, b- 
! 00 


o 3 
45 ^ 

CO r^-J 

„ d 

.2 d 


o 

d 


43 

-H> 


4 


bD 

00 

o 

C/3 


d 

CD 

t-H 

00 

.3 


> 

CM 

t-H 


03 

c5 

& 


43 

-H> 

CD 

bjO 

c5 

a 

CL 

-H> 

O 

-4 

4? 

-r 


4 

‘Bh 

oo 



The total exposure time, shown in units of seconds, is usually the result of 3 or 4 combined (dithered) exposures, with individual times as indicated. 



Date in which the observation took place, and name of the HST observing program that obtained the data. 
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Table 2. Foreground Extinction Correction Values 


Filter 

Band 

k(A) a 

F125LP 

FUV 

8.54 

F275W 

NUV 

6.29 

F336W 

U 

5.07 

F330W 

u 

5.06 

F435W 

B 

4.21 

F550M 

V 

3.05 

F658N 

Ha 

2.54 

F814W 

1 

1.80 

F110W(NIC2) 

J 

1.03 

F110W(WFC3) 

J 

1.00 

F128N 

P P 

0.84 

F160W 

H 

0.58 

F187N 

Pa 

0.46 


a The extinction curve for 
the Milky Way, k(A), ex¬ 
pressed as: F obs (X) = F int (X) 

^g-0.4E(S-y)fe(A)_ The va } ues of 

the extinction curve are from the 
parametrization of Fitzpatrick 
(1999), with total-to-selective 


extinction value Ry=3.1. 








Table 3. Cluster Location and Photometry 
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Table 4. Physical Parameters of the Clusters Outside the Radio Nebula 
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Table 5. Physical Parameters of the Clusters Within the Radio Nebula 
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Table 6 . Intrinsic and Predicted Ha Luminosities 


Cluster 

Log [L(Hg:) intrinsic\ 

Log [L(Hq) predicted] 

(#) 

(erg s' 1 ) 

(erg s _1 ) 

(1) 

( 2 ) 

(3) 


# 1 

qo qo+0.18 

oo.oo_ 0 12 

38-13l°f 

#2 

38.211°;}° 

38.07l°; 3 7 ? 

#3 

37.621°;! 

q 7 77 + 0.72 
' —0.30 

#4 

37.851°;}° 

37 QQ + 0 - 73 

01 ' yy _0.19 

#5 

39.451°;}° 

39 7fi + 0 - 04 
o». i u_o 04 

#6 

<35.36 

37.2611! 

#7 

07 1 rr+ 0.22 

Of. 10 _o 22 

36.81!°;! 

#8 

<37.31 

36-381°;! 

#9 C 

38-621°;}° 

37-461°;! 

# io 

36.81!°;! 

37-561°;! 

# 11 

40.25l°;}| 

40.291°;! 


a The intrinsic Ha luminosity of each cluster, in log 
scale, including 1 a uncertainties, derived from the 
measurements listed in Table [3j corrected for dust at¬ 
tenuation using eq. ( 1 ) and the color excess values 
listed in Tables [4] and [5| For cluster 11, both eq. (1) 
and ( 2 ) are used to remove the effects of dust attenu¬ 
ation. 

b The predicted Ha luminosity, in log scale, from 
the best-fit ages and masses of each cluster, including 
uncertainties (Tables [4] and [5]). Except for clusters 5 
and 11 , the main contributor to the uncertainty in 
the predicted L(Ha) is the uncertainty in the best-fit 
age. For clusters 5 and 11, the main contributor to the 
overall uncertainty in L(Ha) is the uncertainty in the 
best-fit mass, since the ionizing photon rate is fairly 
constant for ages <2.5 Myr. The Ha luminosities are 
related to the ionizing photon rate Q(H°) via: L(Ha) 
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[erg s 1 ] = 1.37x10 12 Q(H 


,-ii 


(Leitherer et al. 


1999 Calzetti 2013). 


c The intrinsic Ha luminosity of cluster 9 is likely 
over predicted, by about an order of magnitude, due 
to the presence of a highly dust-attenuated, younger 
interloper (see discussion in section 6.1). If using the 
color excess derived from the SED fits, rather than the 
line ratios (Table [4]), the intrinsic Ha luminosity of 
cluster 9 decreases to Log[L(Ha)j n f r j ns j c ]=37.23)to:[2, 
in better agreement with the predicted Ha luminosity. 








55 


Table 7. Comparison with the Literature 


Cluster 

Age a 

Age C 97 b 

Age T oi c 

AgeH04 d 

Age dG i 3 e 

(#) 

(Myr) 

(Myr) 

(Myr) 

(Myr) 

(Myr) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


# 1 

C + l 
°-2 

2.5-4.4 


1-5 


#2 

c+l 

2 

2.5-4.4 


1-5 


#3 

5 +4 

°-2 



3-5 


#4 

6±2 





#5 

111 

<3 

n+0.7 

— 0.8 

1-3 


#6 

10±! 

8-12 


10-11 

12-16 

# 7 

10± 4 



6 


#8 

15±t 

10-17 


9-14 

12-16 

#9 

io±i 

30-50 

q+0.9 

°-0.9 

11-14 

6-8 

# io 


50-60 

o+2.5 

8-0.9 

8-15 


# 11 

1±1 






a The ages, and 1 u uncertainty, of the clusters in this 
work from the SED fitting and the EW(Ha) constraints 
combined together. 


’The ages derived by Calzetti, Kinney fc Storchi- 


Bergmann (1994). 


c The ages derived by Tremonti et al. (2001), using UV 
spectroscopy from the HST/STIS instrument. 

d The ages derived by Harris et al. (2004). 


’The ages derived by de Grijs et al. (2013). 
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Fig. la.— A three-color composite of the central 20 ,, xl6" (~300 pc x 250 pc) of 
NGC5253, combining the ACS/SBC/F125LP (bine), WFC3/UVIS/F336W (green), and 
ACS/HRC/F814W (red) bands. The 11 star clusters in this study are identified with ma¬ 
genta circles and numbered. The circle radii are 7.5 pixels (0.19" ~2.9 pc), i.e., 50% larger 
than the baseline photometric apertures used in this work. Clusters 5 and 11 are located 
within the central radio emission region (see Figure [2] and Turner, Beck & Ho 2000). A ruler 
of 1" size (~15.3 pc) is provided at the bottom-left of the figure. 
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Fig. lb.— The same as Figure 
best-fit ages. 


where the clusters’ labels have been replaced by their 


la 
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Fig. 2.— A detail of the region surrounding clusters 5 and 11, using (Left) the same color 
composite as in Figure la and (Right) a three color-composite using the Ha emission line 
(blue), the NICMOS/NIC2/F110W (green), and the Pa emission line (red). The region has 
size 2.4" x 1.8" (~37 pcx27 pc). The two clusters are located at the center of the magenta 
circles, in correspondence of the Ha (cluster 5) and the Pa (cluster 11) emission peaks, 
respectively. The yellow circles identify the nominal positions, relative to the astrometry 
of the LEGUS WFC3 images, of the peak at 1.3 and 2 cm (right-hand-side circle) and the 
secondary peak at 1.3 cm (left-hand-side circle) identified by Turner, Beck & Ho (|2000). The 
absolute astrometric uncertainty of the HST images, 0".l-0".3, is comparable to the radius 
of the magenta circles (0".2). The diameters of the yellow circles, 0".15, are comparable to 
the size of the mm and cm beams (Turner, Beck & Ho 2000 Turner & Beck 2004). The 
orange circle in the left-hand-side panel marks the extent of the 7 mm emitting region, 
1".2 (Turner & Beck 2004), which we term the “radio nebula”. The color-composite in the 
right-hand-side panel highlights: (1) the differential intensity of the Ha and Pa emission 
in correspondence of the two clusters, with Ha being stronger in cluster 5 and Pa being 
stronger in cluster 11; and (2) the color gradient shown by the emission lines in cluster 11. 
The emission map in P/3 is not shown, because of the significantly lower angular resolution 
of the WFC3/IR images relative to the ones shown here. 
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m F275W rn F336W 


m F550M m F814W 


Fig. 3.— Selected color-color plots for the 11 star clusters shown in Figure la The U-B 


(rriF 336 VU~ m ) is shown as a function of both the NUV-U (mF 275 VU~ m F 336 Ub left panel) 
and the V-I ( (m fsso A-r fs i 4 W] right panel). All magnitudes are on the Vega photometric 
scale. Typical photometric uncertainties are shown as thin crosses in the two panels. The 
location of synthetic colors for a range of ages, from 1 Myr to 1 Gyr , is also shown for 
comparison, for both the Padova stellar evolutionary tracks with AGB treatment and the 
Geneva tracks (section 4). Vectors showing the effect on the observed colors of a dust 
correction equivalent to a color excess E(B-V)=0.3 are reported for a starburst attenuation 
curve (black arrow) and for an LMC extinction curve (blue arrow). 
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Fig. 4.— Best fit SED (top-left panel) and 99% confidence histograms for the distribution 
of ages (top-right panel), color excess (bottom-left panel), and mass (bottom-right panel) 
for Cluster 1. In the best-fit SED plot, the blue points with error bars are the observed 
photometry, the magenta points are the synthetic photometry, and the black line is the 
Yggdrasil model, together with dust attenuation and mass normalization, that provides the 
best fit (smallest reduced y 2 value) to the observed photometry. The measured photometry 
in the H-band is not used in the fits, in order to avoid potential contamination by stochastic 
presence of red supergiant stars; we show, however, the predicted photometric value (in 
red) from the best fit model. In the histograms, the continuous lines are for Padova stellar 
models and the dashed lines for Geneva stellar models. The colors indicate: differential LMC 
extinction curve (blue), differential SMC extinction curve (cyan), differential MW extinction 
curve (magenta), and starburst attenuation curve (red). 
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Fig. 5.— Best fit SED plots for (clockwise from top-left) Clusters 2, 3, 4, and 6. The blue 
points with error bars are the observed photometry, the magenta points are the synthetic 
photometry, and the black line is the Yggdrasil model, together with dust attenuation and 
mass normalization, that provides the best fit (smallest reduced y 2 value) to the observed 
photometry. The red point in correspondence of the H-band is the predicted synthetic 
photometry from the best fit model; in all cases, the predicted photometry is within 2 cr of 
the observed photometry, although the H-band data are not used in the fitting. 
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Age (Myr) Age (Myr) 

Fig. 7.— Confidence histograms for the distribution of ages for Cluster 1 (left) and Cluster 10 
(right) for seven bands SED fits (continuous line) and 5 bands fits (dash line). Both a 
differential LMC extinction curve (blue) and the starburst attenuation curve (red) are used 
for Cluster 1, while only the starburst curve is used for Cluster 10. For both cases, the 
starburst curve provides consistent results for both seven and five bands fits, also in the case 
of Cluster 1, which has a better solution with the differential LMC curve. 
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Fig. 8 .— Histograms of the 99% confidence level for the age distributions of Cluster 4 (left- 
hand-side panels) and Cluster 6 (right-hand-side panels), determined using all bands (red 
histograms) and only the five bands from the NUV to the I (black histograms). The top 
panels are for photometry that includes emission lines and fits performed with Yggdrasil 
models. The bottom panels are for photometry from which emission lines have been sub¬ 
tracted, and fits are performed with Starburst99 models. The agreement between all bands 
and five bands fits for the most likely age is higher when the emission lines are subtracted 
from the photometry. 
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Log(X[A]) Log(X[A]) 

Fig. 9.— Best fit SEDs for Clusters 5 (left panel) and 11 (right panel). As in Figure [2j the 
blue points with error bars are the observed photometry, the magenta points are the synthetic 
photometry, and the black line is the Yggdrasil model, together with dust attenuation and 
mass normalization, that provides the best fit (smallest reduced y 2 value) to the observed 
photometry. For Cluster 5, only foreground dust is required to achieve a reduced y 2 ~1. For 
Cluster 11, a combination of both mixed dust, with a total dust column of Ay ~48.7 mag, 
and foreground dust, with Ay ~1.9 mag, are needed to approximate the observed SED. 
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Fig. 10.— Best fit SEDs for Cluster 11, including stellar, nebular, and dust emission, plus 
dust attenuation for the stellar and nebular components. The best fit for the stellar, nebular 
and dust attenuation components are the same as in Figure [9j while the dust emission 
is the combination of two modihed black-bodies, with temperatures T^.i ~ 1,100 K and 
Td ,2 ~440 k, and emissivity e=1.8. The two data points longward of 1.6 //m are from Vanzi 


& Sauvage (2004). Virtually all of the excess in the J and H bands can be accounted for 


with this simple dust emission model. 
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Fig. 11.— The ages, with their 1 a uncertainty, of the 11 star clusters as a function of 
their distance from cluster 11. For each data point, the corresponding cluster number is 
indicated at the bottom of the figure. A systematic trend is observed, with the older clusters 
located further away from the location of cluster 11, i.e., from the region of most active star 
formation in the galaxy. 

















